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ITHE- GBODINS* BESPIBITORY C'ONTECI HOCEL USING 
THE'ONIVAC 1110 EEfiOTE BATCH AHD DEHAMD 

PEOCESSING (General Electric Co’. ), 74 p Onclas 

HC AQ4y MF API CSCL> 06P G3/52 22229 

EROGRAM DESCRIPTION GUIDE ^ ~ 


A. IDENTIFICATION 
Program Name 
Programmer's Name 
Programmer Contact 
Date of Issue 


- Grodins (Respiratory Control Model) 

- Marks , Archer 

- V. J. Marks or G. T. Archer, GE/tSS, Houston 

- 9/6/74 


B. GENERAL DESCRIPTION 

The purpose of this model is to illustrate the transient and steady-state 
responses of the respiratory control system, for variations in volumetric 
fractions of inspired gases and special system parameters . The program 
contains the capability to change workload. 

The program is based on Grodins ' respiratory control model and can be 
envisioned as a feedback control system comprised of a "plant" (the con- 
trolled system) and the regulating component (controlling system). The 
controlled system is partitioned into 3 compartments corresponding to 
lungs, brain, and tissue with a fluid interconnecting patch representing 
the blood. 


C . USAGE AND RESTRICTIONS 


Machine, Operating System, and 

Compiler Required - Univac 1110, EXEC 8, Fortran 

Peripheral Equipment Required - Printer, Card Reader, Graphic Terminal 

Approximate Memory Required - I 867 I 


D. PARTICULAR DESCRIPTION 
See TIR 741-MED-3047 

E. DESCRIPTION OF INPUT 


A Univac 1110 file (GR0DIN.) contains the source and relocatables of all 
the subroutines, the executable program. Batch and Demand data files 
(GR0DIN.GE0IIAT and GR0DIN.GE0DATT) , and Remote Batch run stream (GR0DIN. 
RUNB) , Runs may require a data file and/or the Remote Batch run stream 
file to be modified. Since GR0DIN. is not protected, the user should copy 
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GR0DIN. into another file, and make his modifications to this other file. 
Examples to follow should explain. 

1. Batch - 

The print output will be generated at the 1110 (onsite) . 


a. Onsite Batch 

The fun will be in card deck form, and initiated at the 1110. A 
JSC form 588 a should be submitted with the card deck. 


(l) Job Stream Cards (Start Col. l) 

@RTM, /R DBH0LD , 7007-Q509-C ,DB6-G03432 ,TT, PPP 
(TT=mins .run time, PPP=pages output, MME=user 
name) 

@ASG,AX GR0DIN. 

@C0PI GR0DIW. ,TPE’$. 

@EREE GR0DIR. 

@XQT 

data cards 
©FIN 


^ol.6l 

NAME SKUAB 


(2) Data Cards (See Appendix A for example) . 


Card No. Column Format 


Description 


1 1-4 A4 


Blank if Batch Mode. TTYb if Demand Mode. 


2-49 


Cards 2-49 are initializing cards for 48 
variables (See Appendix A). 


6-20 F15.0 


Variable value . 


26-33 2A4 


Variable symbol (this is not used in the 
program) . 


50 -N 


Cards 50-N are workload cards . 


1-6 f6.2 

10-15 f6.2 


Workload (watts) . 

*Run time (mins) for workload 


N+l 


■5^orkload card with time = 0. 


^Length of computer run is determined when variable on card l6 (maximum length 
of computer run) is reached, or a workload card with 0 run time is read. 


b . Remote Batch 

The run is initiated from a teletype. The run stream and the 
data are in Uni vac 1110 files GR0DIN.RUNB and GR0DIN.GR0DAT. 

The data file is the same format as the Onsite Batch data cards. 
The following example shows copying GR0DIN. into another file 



FILNEW.), modifying the 
run stream into another 

3 

run stream and data files, copying the 
file (FILRUN.), then executing the run 

via a teletype. 


TTY Innut 

Comments 

@ASG,CP FILKEW.,F2 

Catalog a new file named FIUJEW. 

@ASO,CP FILRUN.,F2 

Catalog a new file to put run stream. 
Run stream must be in separate file 


to execute . 

©COPY GR0DIN.,TPF$. 

Copy GR0DIN. into temporary work file. 

©BD,L RUHBjRUNB 

To edit run stream. 

The run stream is the following format: 



^Col. 61 

-©RUN, /R DBHOID , 7 OO 7 -Q 509 -G ,DB6-G03^32 ,TT , PPP 
(TT= mins, run time, PPP = pages output) 
@ASG,AX FIBNEW. 

©C0PY FinTEW. ,TPF$. 

©FREE FILNEW. 

©XQT 

©ADD GR0DAT 

NAME NODECK 


The user might want to change the run time or output pages in this run stream. 


@ED,L GR0DAT,GE0DAT To edit the data file. 

The data file is in the same format as Onsite Batch data cards . The user 
can modify this base data file for his particular run. 

©COPY TPF$. ,FI1K(EW. . Copy temporary work file into FltNEff. . 

©COPY TPF$., FILRUl. To get run stream into FILRIM.. 

©FEIEIE FIIl^TM. 

©START FILRUN.RUTT Start Remote Batch execution, 

NOTE: After above TTY input, FILNEW. contains the latest run ^stream and 

data files. FHiWEW. is not a permanent file, but it might be 
several days before NASA deletes it. The next time the user wants 
to modify these files, he may substitute the first three TTY inputs 
with ©COPV' FIIiNEIW.jTPFI. until FILNEW. is deleted. 
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2. Demand (Time -Share) - 

The run is initiated from a teletype, and the printout will be to the 
teletype. A data file (GR0 DIN.(jR 0LATT) is required to initialize the 
49 input variables. Workloads are then input by the user via respond- 
ing to questions asked by the program. Options are available so all 
workloads can be input before the program executes, or intermediate 
results may be requested before continued input. Rote following 
example’ (Appendix B also contains an example)': 

TTY Input Comments 

@GOPY GRODIR. ,TPF$. Copy GR0DIR. into 

temporary work 
file. 

@ED,L GRODATTjGEODATT To edit Demand 

Data Pile. 


The Demand Data File is in the same format as Batch Data File except 
the first card record contains TTY in first three columns, and the 
49th card record ends the data file (no workload cards). Card record 
16 (maximum length of computer run) and record 4o (time increment for 
printout) are ignored by the program in the Demand Mode. The user 
can modify this data file for his particular run. 

ADD DATA.. . 

@ADD GR0DATT 

The 48 input -variables will then be printed out. 


Execute program. 
Program request , 
User answer. 


IRPUT WORK CARDS... 

W0RK = W0RK L0AD (WATTS ) . . . 

MINS = TIME F0R W0RK L0AD. . . 

PRIRT = TIME IRCEEMENT (MIRS) FOR PRIRT0UT 
fiXSC » » • 

M0RE = IRPUT M0RE BEF0RE EXEC... 

RUN = EXEC. WITH AB0VE, THER CAR IRPUT AGAIR. . . 
ST0P = EXEC. WITH AB0VE, THER ST0P, . . 

BACK = ERASE PREVI0US W0RK REC0RD... 

W0RK MIRS PRINT EXEC (f6 .2,IX,F6 .2,IX,F6 .2, 
0. .2 .1 M0RE IX, A4) 

.00 .20 .10 M0RE 


Program 

outputs 

an 

explanation 
of parameters 
on work 
card records . 


Prog. asks for work record 
User answer. 

Prog, outputs user input. 


W0RK MIRS 
100. 10. 
100.00 


PRIRT EXEC 

.2 RUN 

10.00 .20RUR 


(f6.2,IX,F6.2,IX,F6.2, 

IX,A4) 


Prog. asks for work record 
User answer. 

Prog . output s user input . 


The program executes and outputs the results for the preceding input, 
then (since EXEC was RUR) the program will ask for additional work 
records . 
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W0RK MIUS PRIMT EXEC (F6.2,IX,r6.2,IX,F6.2,IX,A4) . . .Prog.asks for 

work record 

Etc. 

Until the program finishes executing a work record with an EXEC 
parameter = ST0P. 


F. DESCRIPTION OF OUTPUT 

(See Appendix B for Demand example) 

(See Appendix C for Batch example) 

(See Appendix D for normal values) 

TIME (mins) 

ALVEOLAR, ARTERIAL, BRAUJ, TISSUE, V BRAHT, and V TISSUE 

volumetric fractions of C0-, 0 , and K . 

2 ’ 2 2 

ALVEOLAR, ARTERIAL, BRAIN, TISSUE, CSF, V BRAIN, and V TISSUE 
partial pressures of CO^, O^j and N^. (ram Hg) 

ALVEOLAR, BRAIN, TISSUE, and CSF derivatives of the partial pressures 
of CO^, Og, and N^. (ram Hg) 

ARTERIAL, BRAIN, TISSUE, CSF, V BRAIN, and V TISSUE 

hydrogen ion (H+) concentrations (nanomoles) and pH. 

ARTERIAL, V BRAIN, and V TISSUE concentrations of HbO^ (oxyhemoglobin), 
(liters 0 ^ - STED) 

ALVEOLAR RQ 

RQ DIFF 

TRANSPORT TIMES (mins) 

AB = Lung to brain 
VB = Brain to lung 
VT = Tissue to lung 
AT = Lung to tissue 
AC = Lung to carotid body 

VI = Inspired ventilation (liters/min) 

VE = Expired ventilation (liters/min) 

Q = Cardiac output (liters/min) 

FB = Brain blood flow (liters/min) 


DERIVATIVES of Q and FB 
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EESP 5KEQ (breaths /min) 

MINUTE VOLUME (liters/min) 

DEAL SPACE VEtWIIATION (liters/min) 

HEART RATE (beats/min) 

ARTERIAL VENOUS 0^ DIFFERENCE (liters O^/liters blood) 

DEAD SPACE VOLUME (liters) 

CHANGE IN W0RK L0AD (watts) 

METABOLIC RATE CHANGE 0F CO^ AND 0^ C0NSUMPTION IN TISSUE 

At termination of a run, the final values for the 2nd thru 15th 
input variables (see Appendix A) are printed out. If the 38th 
input variable (constant involved in controller equation) 

^ .00001, it is recalculated and printed at the termination 
of the run. 


G. INTERNAL CHECKS AND EXITS 

1. Batch 

The run terminates when time becomes greater than the l6th input 
variable (maximum length of computer run), or when a wor^oad card 
with 0 time is read. 

2. Demand 

Up to 50 worh records can be input before letting the program execute. 
If this limit is reached, the program will execute using an EXEC 
parameter = RUN. The EXEC parameter is checked for validity. The run 
terminates when a work record with an EXEC parameter = STOP is finished 
executing. 


H. INDEPENDENT SUBROUTINES 


See Appendix E for listing of all subroutines . 


I. SYSTEM SUBROUTINES 


No special system subroutines required. Library File TEK. is required to 
supply terminal graphic subroutines. 
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J. COMPLETION OR FIML CHECKOUT DATE 

Appendix A - Data Card Example 
Appendix B - Demand Example 
Appendix C - Batch Example 
Appendix D - Normal Values 
Appendix E - Program listing 



APPENDIX A 

INITIALIZING CARDS FOR 
INPUT VARIABLES 
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Normal Initial 


Card No. 

Variable 

Value 

Symbol 

Description 

1 

2 

ITTY 

c(i) 

.17827 

FA(C02) 

Blank for Batch Mode. 

TTY in 1st 3 cols, if Demand 
Mode. 

Alveolar gas fractions (dry). 

3 

G(2) 

.53459 

FA(02) 

volumetric fraction of gas, 

4 

C(3) 

.28714 

FA(N2) 

dimensionless 

5 

c(4) 

.64121. 

CB(C02) 

Concentration of gas in brain. 

6 

C(5) 

.00116 

CB(02) 

liters (STRD) /liter brain. 

7 

C(6) 

.00105 

CB(N2) 

8 

C(7) 

.61553 

CT(G02) 

Concentration of gas in tissue 

9 

c(8) 

.00147 

CT(02) 

compartment. Liters (STRD)/ 

10 

c(9) 

.00105 

CT(N2) 

liter tissue 

11 

C(10) 

6.00000 

Q 

Cardiac output blood flow, 

12 

c(ll) 

.74913 

QB 

liters /min. Cerebral blood 

13 

C{12) 

48.17427 

PCSF(C02) 

flow, liters/rain. 

Partial pressure of gas in 

l4 

C(13) 

36.69498 

PCSF(02) 

cerebrospinal fluid compart- 

15 

C(l4) 

61.17176 

PCSF(N2) 

ment, mm Hg. 

16 

C(15) 

40.0000 

OMAX 

Maximum length of computer 

17 

C(l6) 

0.0000 

GENT 

run, min. 

A workload card with time=0 
will also end a Batch computer 
run. Demand Mode ignores this 
variable. 

Central sensitivity partition. 

18 

C(17) 

.2000 

SENS PT 
HE 

Weighting of the H+ concentra- 
tion in CSF with that of venous 
blood in the brain. With C(l6) 
=0, zero weight is given to 
venous blood at level of the 
brain and a weight of one is 
given to H+ concentration in 
CSF. 

Blood oxygen capacity, liters 

19 

C(l8) 

.1000 

R1 

(STPD) /liter blood 

Time constants for cardiac 

20 

c(l9) 

.1000 

R2 

output response (El) and cere- 
bral blood flow response (E2) 


for changes in hlood chemical 
composition. 


21 

C(20) 

1 . 1380 . 

GNT SENS 
COF 

Controller sensitivity weight- 
ings; i.e.. 

22 

C(21) 

1.1540 

CRTD HDT 
SCF 
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Card No. 


Normal Initial 

Variable ~ Value Symbol 


V 

1 


^•^38Ccgp(jj+)^ I.I 5 J+O 


(t-l'ao)+ TERM-V , , 
l(N) 


where 

f{ ao 


= Blood transport delay from 


lung to carotid body, 
defined in C(37), and 

TERM = function of 


23 

C(22> 

3.0000 ■ 

KL 

Volumes of lung (alveoli), byain, 

24 

C(23) 

l.OOOO 

KB 

and tissue compartments, liters. 

25 

C(-24) 

39.0000 

KT 


26 

0(25) 

.0500 

MRB (C02) 

Metabolic rates by brain, 

27 

C(26) 

.0500 

MRB (02) 

liters (STED)/min. 

28 

C(27) 

81.9900 

D (C02) 

Diffusion coefficient for gas 

29 

C(28) 

4.3610 

D (02) 

across "blood -brain" , liters(lO)“' 

30 

0(29) 

2. 5240 

D (N2) 

(STH))/min. per mm Hg. 

31 

0(30) 

260.0000 

B 

Barometric pressure, mm Hg. 

32 

0(31) 

.0192 

FI (C02) 

Volumetric fraction of gas (dry 

33 

0(32) 

.7000 

FI (02) 

inspired), dimensionless 

34 

0(33) 

.2808 

FI (N2) 


35 

0(34) 

.1000 

KCSF 

Volume of cerebrospinal fluid, 
liters 

36 

0(35) 

.0000 

T 

Initial time. 

37 

0(36). 

.0078125 

H 

Size of computer time step, min. 

38 

0(37) 

87.5500 

vi(n) 

Constant that is involved in the 
controller equation (See C(2l)). 


Determines the normal level of 
alveolar ventilation so that 
^A(C02) at rest, breathing 
air at sea level. When the con- 
troller sensitivity weightings 
are changed V1(N) should be 
altered accordingly. 


39 


C( 38 ) 5.3900 VI (ss) 


Value used for normal resting 
alveolar ventilation. This is not 
used in the program if V1(N) is 
known. 
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Normal Initial 


Card No. 

Variable 

Value 

Symbol 

Description 

40 

C(39) 

.2500 

FEINT iSIi TIM 

Output printed in these time 
increments (mins). Demand 
Mode ignores this variable. 

4i 

G(4o) 

0.0000 

UNKNOWN 

Importance related to C (39), 
but doesn't seem to be of any 
real significance. 

42 

BC(1) 

.5470 

BHC 03 Blood 

Standard bicarbonate content. 

43 

EG (2) 

.5850 


liters CO (STPD) /liter X, 

44 

BC(3-> 

.5850 


37° where 

45 

BG(4) 

.5850 


X = Blood, brain, tissue, CSF. 

46 

RMT(l) 

.1820 

EWT(C02) 

Metabolic rates by tissue. 

4t 

EMT(2) 

.2150 

EMT(02) 

liters (STH))/min. 

48 

DJ(l) 

.0000 

DJI 

Used in performing Dejours 
experiment (not utilized in 

49 

DJ(2) 

.0000 

DJ2 

present runs). Brief description 


of Devours work relating O 2 and 
CO^ threshold effects is given 

in Grodins' paper. 


In Demand Mode, 

the 49 th card 

record ends the data file . 

If for Batch Mode, cards 50 thru N are workload 

cards . 

Example : 




Card No . Variable 

Columns 

Format 

Description 

50-N W0RK2 

1-6 

f6.2 

Workload (watts). 

DUEA.T 

10-15 

F6.2 

Run time for workload (mins). 


The last workload card should have Run Time = 0. Length of computer run 
(Batch Mode) will he when C(l5) is reached, or a workload card with 0 run- 
time is read. 



APPBHDIX B 

EXAMPLE OF DEMAM) 
INPUT /OUTPOT 
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exoT 

GROOIf^S- RESPIRATORY CONTROL MODEL 


Am DATA . . . 
>@ADD GRODATT 


^.RESPIRATORY CHEMOSTAT — INPUT DATA# 


1 

-1783 

.5346 

.2871 

.6412 

.0012 

b 

.0011 

.6155 

.0815 

.0011 

6.^00 

n 

.7491 

48.1742 

36.6949 

61.1717 

40.0000 

16 

.0000 

.20^ 

.1000 

.1000 

1.1380 

21 

1 . 1540 

3.0000 

1.0000 

39.0000 

.8580 

26 

.0500 

81.9900 

4.3610 

2.5240 260.0000 

31 

.0192 

.7000 

.2808 

.1000 

.0000 

36 

. 0078 

87.5500 

5.3900 

.25^ 

.8000 

'll 

.5470 

5850 

.5850 

.5850 


4^5 

. 1820 

.2150 

.0800 

.0000 



DO YOU WANT GRAPHIC INSTEAD OF TABULAR OUTPUT? (Y/N> 

; 

INPUT WORK i::ARD5. . . 

WORK= WORK' LOADC WATTS > . . . 

MINS- TIME FOR WORK LOAD. . . 

FR I NT= TIME IHCR I MEHTC M I NS >FOR PRINTOUT . . . 

pj/pp 

INPUT MORE BEFORE E}<EC. , . 

RUN = EXEC WITH ABOUE..THEN CAN If'PUT AGAIN. . . 

STOP= E^:EC WITH ABOUE THEN STOP. . . 

BACK= ERASE PREUIOUS WORK RECORD. . . 

WORK MINS PRINT E5<EC ( F6 . 2. IX. F6 . 2. IX. F6 . 2. IX. A4 > . . . 

■••0. 1. .5 STOP 

.00 1.00 .50STOP 


mm-:. 

LOAD CHG.C 

.88WATTS FOR 

l.^INS>AT .0000MINS 


11 ME 

-0O80MINS Pa_U RQ 

8782 RQ 

OIFF 

0027 



ALiEOLAR 

ARTERIAL 

BRAIN 

TISSUE 

CSF 

U BRAIN U 

1 TISSUE 

C02 

,1783 

.5652 

.6412 

.6155 


.6320 

.5998 

02 

5346 

.2021 

.0012 

.^15 


.1352 

.1615 

N2 

2871 

.0010 

.0011 

.0011 


.0811 

.0011 

DKR 

- . 0007 


-.^0 

.^00 

-.0000 




.0018 


.0002 

-.0^ 

-.0108 



[UES 

-.0083 


-.dfm 

".0808 

.0025 



f-'r.D2 

.37.9715 

37.9715 

48.1741 

42.7791 

48.1742 

48.1741 

42.7791 

P02 

113.8677 

113.8677 

36.6162 

46.4015 

36.6949 

36.6162 

46.4015 

PN2 

61 . 1608 

61.1608 

61.1888 

61.1888 

61.1717 

61.1888 

61.1888 

<H+? 


37.6556 

42.3511 

39.0214 

44.0727 

42.9989 

40.0941 

PH 


7 . 4242 

7.3731 

7 . 4087 

7.3558 

7.3665 

7.3969 

HB02 


.1985 




.1340 

.1600 

TRAN- 

i.PORT TIME! 








AB 

UB 

UT 

AT 

AC 

UI 



. 1970 

.1114 

.5912 

.3170 

.1877 

6.0593 



UE 

Q 

FB 

DERIUATIUE8 




c; 9291 

6 . 0000 

.7491 

,0000 

.0003 




PESP FREO 

MIN UOL 

D.S.UENT 

HEART R 

AU02I»=^ 

DSUOL 



12.1685 

7.5753 

1.5811 

66.1070 

43.9640 

.1519 




■^R^GINAL PAGE 1^ 
OF POOR QUALITY 



®J<QT 

G3?ODIHS= RESPIRATORY CONTROL NODEL 


ADD DATA . . . 
>eAOO GRODATT 


^RESPIRATORY CHENOSTAT — IIWT DATA» 


1 

.1783 

.5346 

.2871 

.6412 

. .0912 

6 

.0011 

.6155 

.8015 

.0811 

6.0000 

11 

.7491 

48 . 1742 

36.6949 

61 . 1717 

40.0000 

16 

.0000 

.26^ 

.1000 

.18^ 

1.13M 

21 

1.1540 

3.8008 

1.0000 

39.0000 

.0580 

26 

.0500 

81.99^ 

4.3610 

2.5240 260.0000 

31 

0192 

.7000 

.2888 

.1080 

.0008 

36 

-^78 

87.5500 

5.3900 

.2300 

.0800 

41 

.5470 

.5850 

.5858 

.5850 


45 

.1820 

.2150 

.88^ 

. %JwSfO 



DO YOU WANT GRAPHIC INSTEAD OF TABULAR OUTPUT? <Y/N) 

>Y 


INPUT WORK CARDS. . . 

WORR> WCtRK LOAD<,WATTS). . . 

MIHS= TIME FOR WOf?K LOAD. . . 

PRINT== TIME INCRlMENT<mNS>FOR PRINTOUT... 

e:-:ec 

MuRE= INPUT MOPE BEFORE EXEC. . . 

PUN = EKEr WITH ABOME.THEN CAN INPUT AGAIN... 

STOP= EXEC. WITH ABOUE THEN STOP.. 

BACK> ERASE PREUIOUS WORK RECORD, . . 

WORK MINS PRINT EXEC <F6.2. 1K..F6.2. 1X,.F6.2^1X^A4> 

0 2 . .5 MORE 

.00 2.00 .50MORE 

WORE MTHS PRINT EXEC <F6 .2. 1X.F6.2, IX..F6.2. IX. A4> 

.'■200. 1. .1 MORE 

200. P10 I 00 .10MORE 

WORK MINS PRINT EXEC (F6.2. 1X.F6.2.1X.F6.2.1X.A4) 

. ;:00 . r .5 RUN 

.'•00 00 ?.O0 .50RUN 

TN'PE .SHTFV-OUT C-SO > A^ RETURN— > 



•3 




GRAPHIC aiTPUT(Y.N.S>.TIME INTEffJALS/STARTX/STOP}<AA2/3F5.0>. . 

0 . 10 . 

<A4/8X.F4,0.2F6.0> 


>Y 5. 

P!''nc? Y <srftLF 

PLOKY.K.S) LOC HIGH '*LOW 
:>Y 1. 80, 20. 

pn? Y SCALE 

PLOTC Y.. N.. S > LOC HIGH LCW 
,>Y 1. 140. 70. 

W0RK2 Y SCALE 

PLOKY..H..S> LOC HIGH LOW 
>Y 2. 300. 0. 

c:h+> Y SCALE 

PLOKV..N.S> LOC HIGH LOW 

>Y 2. KKK 


<A4.8X.F4.0.2F6,0> 
» « • 

<A4.8X.F4.0,2F6.0> 

<A4.8X/F4.0.2F6.0) 


>Y 3 . 54 . 42. 

UI Y SC^^E 

PLOTt;Y.M.S> LOC HIGH LOW 
>V 2. 90. 0. 

QB Y SCALE 

PLOT<Y..N.S> LOC HIGH LOW 
>Y 4. 2. 0. 

FREQ Y SCALE 

PLOKY,.H..S> LOC HIGH LOW 
>Y , 5. 50. 0. 

A»;02DF y SCALE 

PLOT'.: Y. N.. S > LOC HIGH LOW 

V 6. 200. 0. 


<A4.8X.F4,0.2F6.0> 

(A4.6X.F4.0.2F6.8) 
« • ^ 

<A4.8X.F4.0.2F6.0) 

<A4.8X.F4.0.2F6,0> 
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GRODINS RESPIRATORY CONTROL MODEL 











ORDINAL PAGE 1^: 
OF POOR QUALITY 


GRAPHIC OUTPUTS Y.. N,. S TIME 
>N 


FINAL UALUES FOR FOLLCWING UARIABLES 


X 

.20746 

r'y 

u. 

,52266 

'7 

26988 

4 

.66543 

5 

00122 

6 

.00100 

I' 

.77849 

3 

.00054 

9 

.00100 

10 

18.96434 

n 

.77073 

12 

48.84239 

11! 

36.28547 

14 

60.42884 
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APPENDIX C 
BATCH EXAMPLE 



APHEHMX C 
BATCH EXAMPLE 


•RESPIRATORY CHeMOSTAT INPUT DATA* 


1 »1783 

.S336 

.2881 

.6413 

• DO 1 2 








A .0011 

.6153 

.0015 

.0012 

6.0000 





0 ^ 



n .7<(96 

48. 1202 

36,6316 

70.6B04 

40.0000 








lA .0000 

.2000 

. 1000 

.1000 

1 . 13BQ 








21 1.154Q 

3.0000 

1 .0000 

39.0000 

• 0500 








2A .OSOO 

81.9900 

4,3610 

2.5240 

260.D0Q0 





a?* 



31 .0192 

.7000 

• 2808 

. 1000 

.0000 








36 .0078 

87.5500 

5.3900 

.2500 

• 0000 








*U .5*(70 

.5850 

.5850 

,5850 









95 .1820 

.2150 

.0000 

.0000 









WORK LOAD CHS. 

( .OOWATTS FOR 

.5QMINS) AT . 

OQOOHINS 







TIME 

.□ooo 







ALV RQ 

.9026 

RQ DIFF 

-.0271 


C02 

02 

N2 

0 E K 

I V A T I V E S 

PC02 

P02 

PN2 

1H + ) 

PH 

HB02 

alveolar 

.1783 

.5336 

.2881 

-.0059 

.0102 -.0043 

37.9779 

1 13.6568 

61.3653 




ARTERIAL 

.5653 

.2021 

.0011 



37,9779 

1 j3t6&6fi 

61.3653 

37.6596 

7,4241 

.1985 

brain 

.6413 

.0012 

.0011 

■ 0006 

-.0024 -.0000 

46.1939 

37.6788 

64.1026 

42.3632 

7,3730 


TISSUE 

.6153 

• 0015 

.0012 

.0001 

-.0003 -.0000 

42.7330 

47,3485 

69.9301 

36,9926 

7,4090 


CSF 




,0090 

.1717 -.9675 

46.1202 

36.6316 

70,6804 

44.0233 

7.3563 


V BRAIN 

.6312 

.1388 

.0011 



48.1939 

37,8788 

64.1026 

43.0743 

7.3658 

. 1374 

V tissue 

.5992 

.1632 

.0012 



42,7330 

47,3485 

69.9301 

4U.D956 

7,3969 

.1617 

transport 

TIhES — 

A8 

va 

VT 

AT ,AC ** 

VI 

ve 

0 

FB 

derivatives 



.1970 

.1114 

.5913 

.3170 .1877 

6.0077 

5,9116 

6,0000 

.7496 

,0000 

-.0025 

RESP FRE8 

10.1710 

MINUTE 

volume 

7.5496 

DEAD Space VENT 

1.5900 

heart rate 

66.1070 



TIME 

.2500 







alv ro 

.9057 

RQ DIFF 

-.0302 


Co2 

02 

N2 

D E R 

1 V A T 1 V E S 

PC02 

P02 

PN2 

«H+) 

PH 

HB02 

alveolar 

,1779 

• 5349 

.2872 

.0016 

.0009 -.0025 

37.898s 

113,9349 

61 . 1664 




arterial 

.5649 

.2022 

.0010 



37.6988 

1 |3,9349 

61.1664 

37.6076 

7,4247 

CO 

> 

• 

BRAIN 

.6413 

.0012 

.0011 

-.0002 

-.0000 -.0000 

48,1942 

38,6339 

63.6564 

42.3634 

7.3730 


TISSUE 

.6153 

.0015 

.0012 

,0000 

-.0000 -.0000 

42.7362 

46,3936 

69.6465 

38,9946 

7.4090 


CSF 




.0087 

-.0002 -.9970 

46.1224 

38.6352 

70.4346 

44,0254 

7.3563 


V BRAIN 

.6321 

.1352 

• 0011 



46,1942 

38,6339 

R 

43,0120 

7.3664 

.1341 

V tissue 

.5996 

.1615 

.0012 



42,7362 

48,3936 

69.6465 

40,0670 

7.3972 

. 1600 

TRANSPORT 

times — 

AB 

V8 

VT 

AT AC •» 

VI 

VE 

U 

FB 

derivatives 



.1971 

.1116 

.5912 

.3169 .1877 

5.9578 

5.8639 

6,0000 

.7471 

.0000 

-.0010 

RESP FREO 

10.1710 

minute 

volume 

7.4971 

DEAD space vent 

1 .5662 

heart RATE 

66. 1070 



work load CHS. 

1 ao.oowatts for 

.2SH1NSI AT • 

5000HINS 







CHANSE IN HETABOLIC RATES HRC02* 

.1892 

HR02> .2150 








time 

.5000 








ALV RO 

.9010 

RQ DIFF 

-.0255 


C02 

02 

NZ 

DERI 

1 V A T I 

V E S 

PC02 

P02 

PN2 

(H+) 

PH 

HB02 

alveolar 

.1781 

• 5352 

.2367 

-.0008 

.0020 

-.0013 

37,9345 

1 |3.9878 

61 .0777 




arterial 

.5650 

• 2022 

.0010 




37.9345 

1 1 3.9876 

61,0777 

37.6313 

7.4245 

. 1986 

brain 

.6412 

• 0012 

.0011 

-.0002 

.0000 

-.OUOQ 

4e. laoo 

36 . 6644 

63.2613 

42.3547 

7.3731 


tissue 

.6 153 

.0015 

.0012 

-.0000 

-.0000 

-.0000 

42.7361 

46,2715 

6 9 , 3 7 0 n 

38.9946 

7, 4090 


CSF 




,0068 

• 0040 

-1 .0180 

48.1244 

36,6357 

70,1826 

44.0271 

7.3563 


V BRAIN 

.6320 

.1353 

.0011 




4&. 1800 

36.6644 

63.2613 

43.0049 

7.3665 

.1342 

V TISSUE 

.5997 

• 1613 

.0012 




42.7361 

46.2715 

89, J7un 

40.0629 

7.3973 

.1598 



APPENDIX C 
BATCH EXAMPLE 


transport times — AB 

.1971 

rESP FREo I0.17ID MINUTE 

ve VT 

.1116 ,59 It 

VOLUME 7«S102 

AT 

,3U9 

DEAD Space 

AC *• 

. 1 B77 
VENT 

VI 

5.9726 

1.5070 

WORK LOAD CH6. 1 lOO.OOWATTS FOR 

.ZSHINS) AT 

•7500MINS 



Change in metabolic rates hrcoz- .*i3A7 

MR02- 

.<(963 



TIME 

.7300 








C02 

02 

N2 

DERI 

V A T I V 

E S 

PC02 

alveolar 

.1485 

.5788 

.2727 

.DI3B 

.0227 

>,0365 

31.6311 

arterial 

*5^83 

,2033 

.0010 




31.6311 

srain 

.6413 

.0010 

.001 I 

-,0'l-24 

-.0003 

- i'oooo 

48.1834 

TISSUE 

,6160 

.0011 

.0012 

,0039 

-.0010 

-.0000 

42.8685 

CSF 




,0068 

-.6030 

-1,0292 

48.1274 

V BRAIN 

,6355 

.1217 

.001 1 




48. 1834 

V TISSUE 

,6071 

.1345 

.0012 




42.8665 

transport 

times — 

AB 

VB 

VT 

AT 

AC •• 

VI 



.1771 

.1200 

.5517 

.2788 

; 1644 

15.8417 

RESP FRE9 

12,3669 

minute 

VOLUME 

18.9704 

DEAD space VENT 

2,6828 

time 

1 .0000 








C02 

02 

N2 

DERI 

V A T I V 

E 5 

PC02 

alveolar 

. 1S36 

,5774 

.2690 

.01-08 

— . 0 2 0 2 

.0094 

32.7114 

ARTER I AL 

.5350 

.2032 

.001 0 




32,7114 

BRAIN 

,6379 

.0010 

.001 1 

».01 10 

.0000 

-.0000 

47.4443 

TISSUE 

,6161 

,00DB 

.0012 

,0006 

-.0005 

-.0000 

42.8917 

CSF 




>.0820 

■ -.6229 

- 1 , QS04 

48, 1 176 

V BRAIN 

.6323 

.1212 

.0011 




47.4443 

V TISSUE 

.6146 

. 1052 

.0012 




42,8917 

TRANSPORT 

times — 

AB 

VB 

VT 

AT 

AC *o 

VI 



,1472 

,1195 

.4364 

.2184 

. 1345 

17.9397 

rESP FReq 

15.2405 

MINUTE 

volume 

2 1 .4438 

DEAD Space vent 

3. 1 503 


VE 9 

FB 

derivatives 

5,8738 6,0£j0n 

.7480 

,0000 ,0005 

heart rate 

66. 1070 




ALV Kl) 

1 .6623 

R4 DIFF 

-.7868 

P02 

PN2 

(H+ ) 

PH 

HB02 

173,2812 

58,u877 




123,2812 

58,0877 

33,3900 

7.4764 

.1994 

32,1767 

62,9292 

42, -3668 

7.3731 


35.0427 

69,0534 

39,0773 

7.4081 


3i,5570 

69,9266 

44,0299 

7,3563 


32,1767 

62,9292 

42.7587 

7.3690 

.1207 

36,0427 

69,0534 

39,6795 

7.4014 

.1334 

VE 

4 

FB 

derivatives 

l5,7335 

7,3924 

.6117 

1.5630 

-.1929 

heart rate 

78,4256 




ALV RO 

1.3762 

RS OlFF 

-.5007 

P02 

PN2 

( H+ > 

PH 

HB02 

122,9327 

57,3059 




172.9827 

57,3059 

34. 1311 

7.4668 

.1993 

3l .8706 

62.3184 

41.9056 

7.3777 


76,4801 

68,5792 

39,0919 

7,4079 


36,3957 

69,6638 

44.0210 

7.3563 


3 t . 67 o6 

62.3)64 

42.2958 

7,3737 

. Uo2 

?6.480l 

68.5792 

39,1903 

7. 4068 

.1044 

VE 

4 

FU 

derivatives 

18 , 6461 

9,2192 

.6102 

• 2.7336 

.0437 

heart RATE’ 

94*5195 





APPENDIX C 
BATCH EXAMPLE 


FINAL values for following VARIABLES. 


1 

. 15357 

z 

.57738 

3 

.2690*f 

N 

.63787 

5 

>00101 

6 

.00107 

7 

.61609 

S 

iOooei 

9 

.00118 

lo 

9.21921 

11 

>61019 

12 

98. 1 1769 

13 

36.39575 

I*! 

69.66389 



APPENDIX D 

OUTPUT VAEIABLES 
NORMAL VALUES 



I 


Normal Variables Values are for conditions: 
B = 260 ram Hg 

^I(C02) = •‘^^92 
^1(02) 


Physiological 

Variable 

Computer 

Variable 

Normal 

Value 

Phys iologi cal 
Variable 

Computer 

Variable 

Normal 

Value 

^ACcOg) 

c(i) 

.1783 

^CSP(C02) 

C( 12 ) 

48.1743 

^A(02) 

C( 2 ) 

.5346 

PCSP(02) 

C( 13 ) 

36.6950 

^A(N2) 

• c(3) 

.2871 

^CSP(N2) 

C(l 4 ) 

61.1718 

P / 

a(C02) 

F( 7 ) 

37.9705 

^a(C02) 

CC( 1 ) 

.5652 


P(l) 

113.8681 

‘^a(02) 

F( 9 ) 

.2021 

^a(N2) 

PAN 2 

61.1614 

^a(N2) . 

_P( 10 ) 

.0010 

^B(C02) 

CPB 

48.1745 

^B^COg) 

c(4) 

.6412 

"■bCOj) 

P( 17 ) 

36.6954 

^ 3 ( 0 ^) 

C( 5 ) 

.0012 


F{ 18 ) 

61.1616 

^B(N2) 

C(6) 

.0010 

^tCco^) 

CPT 

42.7786 

‘^1(002)' 

0 ( 7 ) 

.6155 

^1(02) 

PT 02 

46.2605 

^T( 02 ) 

0 ( 8 ) 

.0015 

^TdTg) 

PTN 2 

61.1673 

^TCNg) 

0 ( 9 ) 

.0010 

^a(Hb 02 ) 

CHB(l) 

.1985 

"^aB 

AB 

.1970 

^VT(Hb 02 ) 

CHB( 3 ) 

.1342 

'^aT 

AT 

.3170 

^(HbOg) 

CHB( 2 ) 

.1597 

^vB 

VB 

.1114 


2 


Physiological 

Variable 

Computer 

Variable 

Normal 

Value 

Physiological 

Variable 

Computer 

Variable 

Normal 

Value 


CH( 1 ) 

37 . 6^49 

TvT 

VT 

.5912 


CH( 2 ) 

42.3514 

'Tao 

AC 

.1877 


CH( 3 ) 

38.0212 

VI 

VI 

6.0586 

^gsf(h^) 

CH{ 4 ) 

. 44.0728 

VE 

VB 

5.9248 

^VB(H+) 

HVB 

43.0032 

Q 

0 ( 10 ) 

6.0000 

^vtCh"^) 

HVT 

40.0892 

FB 

C(ll) 

. 1 ^ 9 i 

PHa 

CPH(l) - 

7.4242 

EESP FREQ 

FREQ 

12.1637 

P^csF 

PHCSF 

7.3558 

MimJTE VOL 

TVKT 

7.5723 

^“VB 

PHVB 

7.3665 

DEAD SPACE VETO) 

DEADVT 

1.5806 

p^Vr 

PHVT 

, 7.3970 

HEART RATE 

HRATE 

66.1070 

^Seaiw 

CPH( 2 ) 

7.3731 

ALVEOLAR RQ 

RQ 

.8754 

^^ISSUE 

CPH( 3 ) 

7.4087 

RQ DIFF 

QF( 5 ) 

.0001 

AV 02 DF 

AV 02 DF 

44.1667 

WORK LOAD 

W 0 RK 2 

0 . 

DSVOL 

DSVOL 

.1518 

MRTCOg 

EMT(l) 

.1820 




MRTOg 

RMT( 2 ) 

.2150 


APPENDIX E 
PROGRftM LISTING 





dimension CISD), XM‘»0»2)» SVfl8»50), VTRANC18), RK(H»H), 
i~ TCrr‘liS). DCnM)! A16J| D«I5 >i F(2C!)i VOLIIGJi RMT(2)f' 

2 RCM*i QFl6)t TAU I 5 > I ' CC t 3 ) , ChB(3). CHM>i tPH(3)» 

g- RQThT 

C I 80 ) 


alveolar VOL GAS functions 
1 F A ! C02 J 


F A ( 02 ) 

F A { N2 ) 


GAS concentrations IN BRAIN. 


*1 CB(C021 

5 CB (02 I 


6 CB ( N2 ) 


GAS concentrations IN TISSUE. 
7 CT{C02> 


CT (02 ) 
CT ( N2 ) 


cardiac OUTPUT. 



**I C controller equation SENSITIVITY WEIGHTINGS* 














a^OD > in j: u|N)«^ o<«a9 xi c^ lh xui m^ o 


DN2 


sa c 

S 9 c ba'rometr i c pr’essukeT. 

AO c 3n H 

Af ■ c vol.fract f on"of inspired gas. 

A2 C 31 F I ( C02 » 



c 

32 

F I (02 I 






C 

33 

F 1 ( N2 1 






c 

c V n l.df cs f » 


C 3*4 KCSF 

c initial tihe 


C 

C 

35 T 

COMPUlFR TIMF STFP. 


C 

C 

3A H 

controller Equation constant ( ma i nt a i ns resting 

PA(C02» APPROX. qn}. ■ 

C 

C 

37 9 ( {N 1 

VALUF for resting alveolar ventilation. 


c 

3a 

V 1 (SS » 


c 

output 

PRINT increments (ALSO PRIMS AT .SHIN. 

INCR IHENTs ) . 

c 

c 

39 

PR I NT- ALL time 


c 

sv ( 1 a ,50 ) 


c 

arterial gas concentrations at lung exit. 


c 

1 

CA { C02 J 


c 

2 

C A ( 02 ) 


c 

c 

3 

C A 1 N 2 ) 


c 

VENOUS 

fiAS CONCENTRATIONS AT BRAIN EXIT, 


c 

9 

CVB ( C02 ) 


c 

5 

c vg 1 02 1 



C h CV9( N2I ^ 

C VENOUS GAS CONCENTRATIONS AT TISSUE E X.'I T . 
~ 91 C 7 c'vf(C02J 

92 a CVT(021 

93 C 9 CVTIN2J 

VS C 

"9S ~ c “cardiac output. 

9A C 10 

'97 C CERFPRAL 0LOOO) FLOW, 

9a C I 1 QB 

99 C TISSUE PILQOD FLOW. 

100 C 12 UT 

"l^T ’ ■ C ARTERIAL H + “ CONCENT RAT I ON'. ■ 

102 C 13 C A ( ) 

TtT3 C — AlTTETniSL“T3'2“TF'irSTW; 

1 OH C 19 P A 1 02 1 

17^5 C ^ 

lOA C 15 

1117 C TOTAL GAS CONCENTRAT I ONS" aT "B H A I N EXIT. 

108 C lA CV0(Cn2) + CVei02> + CV«(N2) 

"TUT C TOTAL GAS COli C/rN T K A T IONS AT ~T 1 5 S U E E X lT7“ 

lin C 17 CVTIC02) + CVT1021 + CVT(N2) 

Til c time. 

112 c la t 









rni 

cr 

VTRANl 181 



US 

C 

arterial gas concentrations at brain entrance* 



t 1 6 

C 

1 

CABlcnZl ■ CA1C02MT - TA9» 



U7 

C 

2 

CA8(Q2) » CAIQ2) (7 - TAB) 

' 


118 

C 

3 

C AB I n2 ) « CA(N2)(T -TAB) 



.119 

C 





— r?n 

(T~ 


brain gas CONCENTRATION" T LUNG ENt^ANCE. 



121 

C 

8 

CVB ( C02 ) ( 7 - T VB ) 



1 z z- 

C 

5 

rVBTo7TTT~=~TT¥ ) 



1 2.1 

C 

6 

CVB ( N2 ) ( T - .TVb ) 



1 2R 

c 





1 2S 

c 

VENOUS 

TISSUE GAS concentration AT LUNG ENTRANCE* 



12« 

c 

7 

CVTICnZ) (T - TVT ) 



127 

c 

8 

CVTI02 MT - TVT ) 



1 28 

c 

9 

CVT ( N2 ) ( T - T VT ) 



129 

,c 





130 

c 

arterial gas concentrations at tissue entrance. 



131 

c 

10 

CATIC02) *> CA(C02)(T - TAT^ 



132 

c 

1 J 

CAT < 0? ) o t A ( 02 1 ( T - T AT ) 



133 

c 

12 

CAT(N2J a CAIN2) {T - TAT) 



S 38 

c 





135 

c 

arterial h+. concent rat ion at carotid B0DIES»SITE* 



13A 

c 

13 

C AO (H+ ) » C A lh+ ) n - T AO ) 



1 37 

c 

arterial 02 TENSION AT CAROTfO BODlES’SITE. 





PAO{02) » rA(02)(T - lAn) 

arterial h+ concentration at hrain entrance. 


G E 


16 tCVB(COZ) ♦ CVB(02) + CVB(N2))IT - TVB) 

total gas' concentration from tissue at lung entrance. 






















- 


173 

C 

7 

P A ( C02 ) 

17M 

c 

8 

P A ( 02 > 

175 

C 

9 

CA t 02 J 

176 

c 


C A ( N2 ) 

177 

c 

1 1 

CA(C02) + CAt02) * CA(N2) 

178 

c 

12 

CVH (02 ) 

179 

r 

1 3 

CVT (02 ) 

180 

c 



1 6 1 

c 

PRODUCT 

OF DIFFUSION COEpFS.ANO GAS DIFFERENTIALS ACROSS BLOCD-BRAIK 

182 

c 


• 

1 83 

c 

19 

DC02 IPB(C02) ■« PCSF(C02)) 

1 89 

c 

15 

D02 (Pa (02 ) - PCSF (021) 

1 85 

c 

1 6 

DN2 ( PB ( N2 J - PCSF ( N2 ) ) 

I 86 

c 



187 

c 

1 7 

PB ( 02 ) 

188 

c 

18 

PB ( N2 ’ 

1 89 


DIMENSION XNB(9,2)» Dj(9), IDJ(2| 

190 


CnMMON/27 r, XN» SV, vTRAN, RK, SC» DC, A, C. F. UOL, RMT. Rr. OF. 

191 


1 

TAU, CC, CHB, CH, CPH, D«, VE, VI, CPB, CPT» CACK, X, DT, 

1 92 


2 

IRK, LOC, ITERX. INDEX. I, J. M, N 

193 


common /R/ XDS,XMH,CXT,WORK,OUMl,OUM2,DUM3,ttORK2,RMTB,RMTB2»TlMEOF 

! 99 


! iRML I N , I n Y 

195 

c 

ITTY»FL<iFORTTYMODE. 

196 

c 

Oa OUTPUT TO PRINTER (BATCH fionEJ.' 

197 

c 

*TTY 

*= TTY. I/O AND 1ST TIME TO SUBROUTINE RC12* 

198 

c 

1 ■ 

TTY I/O and not 1ST TIME TO RC12. 

1 99 


data 

ITTTY/’TTY •/ 

200 

c 

data 

FOR INITIAL CONDITIONS 


201 WRITE UjS) 

202 5 format </» G RODINS; RESPIRAT OR Y C U ia.R OL MO D.EU /.JL, 


203 300 CONTINUE 

20R WR 1 TE { 6 , H83 ) 

2TT5 *»83 FORMAT (• OADU DATAf*#*) " 

20A C read indication DC BATCH OR TTY MO^E. 

lOT REXdT 5T^801 fTn 

208 Rao FORMATCaR) 

20R IF (I TTY .NE. IT TTY) nTY'“^“Q 

2*0 WRITE<A»‘?0) 

2 1 1 90 FORHAT ( 1 H 1 i 1 X » 3 7 H • R ES P I R A T 0 R Y CHEMO^AT w INPUT DA TA*/ ) 

212 C DATA FOR INITIAL CONDITIONS 

2T3 bo 10 I “ 1 I HO 

21R C IlOA HAS PROBLEM WITH END«» > SO THIS 15NT USED TO 

2 T 5 C determine end OF KUNINOO CAPAHILUT T 0 “ STAR T ANOT her 

216 C model run in SAME COMPUTER RUN). 

7T7 Re a d ( 5 rr^FrEN b » 3 o i j c ( r ) . Txrrrn u i , 0 = 1 , 2 )' ‘ ' ' 

21 s 10 CONTINUE 

C ESTABLISH COMPUTER STEP INDEPENDENT OF INPUT DATA* 

C ( 36 )»«76 1 25E-2 

190 FORMAT (5X»F 15.0, 5X ,2AM ) 

00 20 I ° 1 >*) 

IPRO « f + RO 

read t5,19Q) 8C<I1, (XNRtl.J), J ■ 1,2) 

20 continue 

no 30 I “ 1 « 2 


read (5,190) RMT(l), (XN8(I,J),,J ■ 1,2) 


















228 



1P9H « I ♦ 99 

229 


30 

CONTINUF 

230 



00 90 1 = 1T2 

231 



RFAD (5,190) Og(IJ, (XNR(I,J)i.J * Ji2) 

— ?X2 



IP90 a T ”"91i 

233 


9n 

CONTINUF 

239 

C 



235 

c 

nUTPUT INPUT DATA* 

234 



J * 1 

237 



DO 75 I *> 1 » B 

238 



J X a J 9 

239 



VHR ITE t 6 i92 ) J. IC( 12) . I2 = J,0X) 

2‘»0 


92 

format ( * ' 1 12, 2X ,5Cf 9,9) > 

291 



J ■ 0 + S 

292 


75 

continue 

293 



HRITE(6,92) J , « BC( J ) , 1* 1 ,9 ) 

299 



0 " 9 5 

295 



»RITE(6,92) J,HMT!l),RMT(2),DJ(I)iDJ(2) 

296 

c 



297 



COMMON/ PLTHUF/NBUF,ABUF( 181 ),Y8UF( 181, 8), NA(B),KSTOPp,TRAX,KPLT 

298 



DATA KT/IHY/ 

299 



TMA XoC (15) 

2^0“ 



I TE (6,921 

251 


92 

FORPAT(/*ODO YOU WANT C-RAPhlC INSTEAD OF TABULAR OUlpUT? (Y/N)*) 

252 



read (5,99) KY.Y 

253 


99 

FORMAT ( 1 A 1 ) 

259 



IF (KYY,EQ.KY) KPLT»1 

255 

c 

I F 

TTY I/O MAX«TIMt VilLL COME FROM ftORK CARO, 

256 



IFtITTY ,NE. 0) C(lb) « 9999999999. 

257 

c 



258 

c 

FI <CQ2) 

259 



OUM loC ( 31 ) 

260 

c 

F I (02) 

2&l 



0UM2=C ( 32 ) 

2*2 

c 

FI (N2 ) 

2*3 



DUM3»C(33) 

2*9 



ttORKaO . 

2*5 



»>ORK2*0 . 

2*6 

c 

r'tTABOLIC RATF OF 02 CONSlHPTiON I TISSU^i 

2*7 



RMTBaRMT(2) 

~T^V— 



RMTB2®R^^T C 2 i 

2*9 

'c 



27T1 



7 THEOF »n. 

271 



XOSaO. 

i 7 2 




273 



MMMaO 








4 - 


CCiZ) s C {*» > 

TISSUE. CONCENTkAT 1 Oh OF C02. 


C C I -D ) = C ( 7 ) 
brain C02 TENSION. 


CPB n BO.O 
TISSUE C02 TENSION. 


PAGE 1% 

OF POOR oliaiiTv 


SETS VARIOUS 
TMA X C 


Cl 15) * C( 15) 
PRINT ALL TIME. 


VARIOUS constants AND AGGREGATES 


CONSTANTS 


.000 1 


00 200 I » 27,29 
C t I ) a cm * l »E-7 


200 CONTINUE 
!02 continue 


IRK 3 t 

M a I H 

N a 5 

inu ( I ) = 0 


SOLUBILITY COEFFICIENTS. 


AIM* 


AIR)* 

All) 


A I 2 ) 
AO) 


AM) 
A 15) 


( alpha ) C02 * 


( ALPHA )C02 I 
» 0.51 
a C 70 2*1 
=• 0.0 13 
■» 0.5 1 
“ 0.02*1 


A t 2 I » 


A I 5 ) • 


I alpha ) 02 I 


{ alpha ) 02 » 


A O ) » 


A I 6 )« 


( ALPHA) N2 » 


( ALPHA)N2 


1 c 

A ( 6 ) a O.Q O 

ATM/M MH 6 COMVEhSION FACTOR. 


SK a 0 , 

.00132 

1 c 

carbonic ACID DISSOCIATION CONSTANT,. 


C AOK a 

795.0 

1 c 

VOL ( } ) "VOL ( 1 0 ) = VOLUMES USED IN CALCULATION OF VARIABLE TIME DELAYS. 


VOL ( 1 ) 

= 0.015 


VOL ( 2 1 

a 1.062 

1 

VOU { 3 > 

® 0 « 1 

1 

VOL ( B ) 

* 0 « 0 6 

\ 

VOL(S) 

= 0.188 

I 

VOL ( 6 ) 

a 2 .99 


VOL t a ) • 
VOL (9 )“ 
VOL ( 10 )' 


1.062 


0 .008 
1*062 


{METABOLIC RATE OF C02 


same for 02 












in -c 



D«9) n 863 »Q/ 0 < 1 ) 

C FACTOR USED IN ESTA 8 H 5 H 1 N 6 CB(C 02 ), 


D) «• 0.6 

c manipulation of computer time step* 


DMM) a c( 36 )* 2.0 

Dtl 5 ) a 0 ( 1 *»> -> .Ql*C( 36 ) 


! C 

1 CALL 

RC3 



1 CALL 

• CALL 

RC9 

RC5 

(CPB, F(9) , CM) , 8C<2» ) 


• CALL 

RC21 

tCH8(2), F(3), FlN)| C(R), CH ( 2 ) . 

CPH ( 2 ) ) 

' CALL 

RC 1 9 

(CP8* CH8I2)* CC(2), acin* F(R)) 


! CALL 

RC5 

(CPT. F <6) . C(7 ) , ac<3) ) 


' CALL 

RC21 

(CH813)* f<5), PU), C(7)» CH<3)i 

CPh t 3 ) ) 

) CALL 

RC !9 

tCPTi CH6(3)i CC13), bCU), F(6)l 


CALL 

RC20 



! CALL 

RC7 



1 call 

RC9 


• 


CAUL RC 9 
CALL RCin 


RC 1 1 
CALL RC 12 


CO TO 60 
50 CALL RClS 


CALL RCI 6 
60 CALL RCI 3 


CAUL RC 12 


80 IF (KPLT.LT*!) 00 TO 76 
KSTOPP * 1 


76 AR1TE<6,78) 

FORMAT t *l~riNAL VALUES FOR FOLLOVli N'O VARIABL E S , ■» ) 
IF (C( 37 ) .<iT. l.QE- 5 ) GO TO' 250 


TfcHM B Q ,0 
IF <VTRAN{l<n - lOH.Q) 230* 2ND, 2 


OH.Q - VTRANI 1 ‘n ) ••*) ,9 J 
2H0 C(37) » C(20»#(C( 16J*VTRANt 15) ♦ (1*0 « C ( 1 6 ) ) *CH ( R ) ) 


I « 37 


250 DO 260 T » 1 * IH 










BR I TE t 6 »83Q >' 

830 format ( * 0HQRHAL TERRATlON*) 

301 continue 

5T08 

C FORMAT '( IH I R8X37h«ftESP IRATORY CHEmOSTAT -- INPUT DATA*///) 

C 92 FORMAT ( R 2 1 3 i 1 OXF I 0 , R , I Q X 2 A& ) 

C 190 FORMAT ( 5 X F I 5 . 0 . 5 A 2 A 6 ) 

192 FQRMAT(* » t 13 >2* »F 1S»5 »2X ,2AR ) 

198 format ( 1H 1 ) \ 

END 



1 ' StjBRQIj-T tNE KC3 

2 triMFNsroir'c r MDTT“irN t *♦ o » 2Ti s v ( i a » 5 qT‘,~vT^a‘ n TTF rr"^ kT iT tTiT, 

3 1 SC(IM,5). DC(l*t)i A(A», D<15)t F(20)» VoL<103» Rf'T(2). 

" 2 “ CTF-,— ITF'ri ) ,~TAU(5) . CC(3), CHB(3), ,“Tp h'( 3 ) ' 

S 3 DQ < 4 ) 

COMPON/2/ C» XN, SW, V7RAN* RK, SC* DC» A,“D, F» VOL, R M f7~B'c ^"QfT 

__J T AU I CC> CHBt CH', CPH, oa, _V JIP B ♦ CPJ_»^CACK_, X, OT_,_ 

‘ “ “2 ' ■ 1 R K r L 0 “C’i HERXi '! NUEX , ’i V’ 0 » B . ' 

C69 A 9 _ ^0 R M A t ( 1 H 7 HS<JB RC3) 

c 'seTs~'i i beViTe^WlnT expressions ' 

C TISSUE blood FLQiA* 

■ (5F(1) *> c ( ill ) ~ 

C ARTER I AL 0 2 T.E MSTON, ^ _ 

i F{'n a”o( r» «c ( ?. ) 

C arterial C02 CONCENTRATION. 



16 F(2)atM5)*C(l) 

17 C BRAIN 02 C C NC F- 14 T R AT J ON / < C ON V • F A C T OR • S CLUB I L I T Y C0EFF*FCB 02) 

TR F(3) * C(B)/D{3J ~ ^ 

19 C ICONv»FACTnR*SOLUBlLITY COEFP.FOR L02) • BRAIN C02 TENiilCNt 

TB ■ Ff 4)“:^DT2'’*CP0 ■ ~ 

21 C TISROE 02 CONCENTRATION / { C 0 N V . F A C T 0 R • S OL U R I L I T T CQEFF.FOR 02) 

2 2 " ■ F(5) » (f(0V7rT(3) 

23 C fCONV.FACTOR«SOLUBlLlTY COEFF»FOR C02) • TISSUE C02 TENSION, 

24 F (6) = n{2)*CPT 

25 C arterial C02 TENSION. 

'26 ‘ F 17) ="b( 1 '«Ct'5 )' 

27 C arterial 02 TENSION. 

28' • - -- F(8T »■ 0« 1 I •’C<"2 ) 

29 return 

_ __ 






— - 10 


OF POOR 


QUALITY 







PAGE \P:. 
M POOR qual:t' 


tJ&-Q03‘»32*Tpl^S»RC5 

I SUBROUTINE RC5 (CP» FB» CCH , 8HC} 

dipensToT) TT^^r ^nmo, 2I, sviieiso). vtran( lei . rk < m *<) > , 

! SC(IH,5)| OCU«J» A(6)| D(15)i E(20>i VoLMOf RKT(2), 

2 n Cl R J I Wf VATV tau ( 5 > f ' CC < 3T , "CHBT3 ) • C H < R) , CPH ( 3 ) » 

3 OQ { R ) 

COMMON/2/ Cl XN* 5V, VTRANt RK , SC » DC i A» D, Fi VOC, RPJ. BC, GFi 

1 TAU. CC. CHB, Ch, CPH, 0^, VE, VI, CP0 , CPT* CADK, X, OT, 

2 "irk; LOX.'^^ERX, I’nETEX, 'l, J, N, N ■ ■ ■■■ 

C iterates for brain and tissue PC02 
X5X69 "■ H 7HSUB RrS) 

510 X B {CCB - FB)/tO.Ol«CP) 

X « RCF1(X> ^ ^ 

C SEE EQUATION R • 1 , X « PB(C02» , 

CCB + D110)*IX - 0,lR))/0(2) 

C CPaPFI(C02). 

■ C ALL^Rr^rcF") ■ 

CP = CP + (X . CP 1/10,0 

18 C cerebral blood FLOW, 

19 FB a 0 ( 2 ) *CP 

20 C3O00 FORMATllH »RHCP*= ,El6,6,RHFB3 E 1 « • 6 t 5 HC’C0» E j 6 , 6 , 5 MBHC* E16»6) 

21 IF IITERX) 520 , 5 1.0, 520 

'2T' ~ “520 return ^ 

23 end 

. 


PRTtS RCA 


2 

3 

R“" 

5 

— s — 

7 

-gp-- 

9 

“TU 

1 1 
— T 2 

1 3 

-yif- - 

15 

TS~' 

1 7 





B6-603‘<32*TPFS*RC4 

1 subro u tine NC6 (Y) 

OIMENS ! ON “C I MO ) .. XN ( MO ,2 ) , SV(18»50I| VTRAN(18)* RKHM»M), 

1 ■SC(IM»5>, DClIMlt A<61t p(i5l> F(20)» V0 lOQ)» RNT(21» 

0C(M), (}FU)i TAU(5»i CCU), CHB<3)* CH(M», CPH(3 >i 

D*a ['M ) 

COMMQN/Z/ C» XN, 5V, VTRANo RK, SC» 0C» a, D, F, VOL, RMT, BC» QF, 
J TAm CCi CHBi Ch, CPHi DSi VE,, VI, CPR, CPT» CADK, X, DT . 

2 IRK* LOC, ITERX, INDEX, 1, J, H, N 

CHECKS CONVERGENCE OF ItERATlVE PROCEDURES 

CM ; X3CA(C02), T«CC(1) • 

C5 5 XaPB I C02 ) I T«CP • 

RC19 : X»CVa<CD2l, Y-CVC , 

:6969 FORHAT i >H 7HSUB RC&) 

I TERX a' 0 ‘ 

DIFF a A8S < ( X - T > /Y ) 


U IF IDIPF - l.OE-51 620, 620, 630 

17 620 ITERX = I 








1 subroutine RC7 

Z DIMENSION C^HQlt XM‘JO/2'), 5V(1.8i5D>i VIRANdS), RK(lHiH), 

3 I SCCntSJi Dcm), A<A>i DJISJi F(2Q), VOLUOJ, RKT(2)t 

^ 2 8C('l)i QF(6). TAU(5)» CC<3 )i CHB(3)f CPH(3), ~ 

5 3 DQ{‘U 

S COMMON/Z/ C, XNi SV , V7RA-;i RK , SC> DCi A, D, F* VOL, RMT. BC» <}Fi 

7 ! TaUi CC» CH8> CH> CPH* DQi V£ » Vli CP8» CRT i C A DKt Xt DT t 

8 2 IRKiLOCiITERXi INDEX* liJiHiN’ 

9 ^ CQMMQN /R / XOSiXhh*CAT,lftORK«DUMl*DUM2,OUM3,VtORK2*RMTfi*RKTfl2*TlwEQF 

10 I *rhli’n’ 

1 1 CA9A9 FORMAT ( Ih ZHSDO RC? ) 

12 C FILLS SV array niTH INITIAL CONDITIONS 

13 CALL RCI& 

P4 IF t XDS ,GT . XMhT (jOT02 

15 00 725 I a 1»1V 

16 DO 720 J B 2i50 

17 SV n I J ) » SV ( I 1 1) 

Te 720 CONTINUE 

19 725 CONTINUE 

2D 2 continue 

Z-1 DO 730 U a 2*50 

T2 SV I 1 8 , J ) a SV I 1 8“,“J - 1 ) . D ( I 9 ) 

23 730 CONTINUE 

C3000 FORMATMH .12H18SV S D I J M) * 6 ( 3 X * E t 6 ♦ 6 ) / I‘H . 6 { 3 X , E 1 6 ♦ 6 ) / I H *7(3X,E1 
C C6*6 ) ) 

return 

END 




B6-G03*132*TPFS.RC8 


1 


subroutine RC8 

2 


OIMFNSION C1*10»i XN<*10*2>, SV(IB,50>, VT.RAN(18), RKll*!,*!), 

3 


1 

SC<l*ti5J» DC < 1 H J ». A 1 4 ) » DllS), FI20I* VOLUO)# RMT(2), 

*1 


2 

BC(«1), (3F14J, TAU(5», CC(3), CHBI3)* CH(*11, CPH(31* 

5 


3 

Ot( t 8 ) 

6 


COMMON/Z/ c, XN, sv, VTRAN* RK, SC* DC* A, D, F> VOU, RMT-, BC, 61F» 

7 


! 

LlMgSKirmi— 

8 


2 

IRK* i.nc, iterx* index* 1* J* M. N 

9 

C 

calculates TRANSPORl TIMES 

1 0 

C 

equations 8, 10 THRU 8.1*1 . 

1 1 

C6969 


format 1 1 H 7HSUB Rc8 ) 

12 


DO 

870 I ■ 1,5 

J 3 


DT 

■ C«35) - SV(ia,l) 

1 *1 


NO 

o 1 

15 


GO 

TO 1810*812,81*1,816,810' . 1 

t 6 

8 1 0 

NC 

» M 

1 7 


NR 

«* 1 n 

18 


Go 

TO 820 

1 9 

8 1 2 

NC 

■ 1 0 

20 


NB 

■ 1 1 

21 


GO 

TO 820 

22 

8 1 *1 

NC 

a 1 0 

23 


NB 

B 1 2 

2*1 


60 

TO 820 

25 

816 

NC 

■ 12 

26 


N0 

■ in 

27 


QA 

« qf U ) 

28 


GO 

To 822 

29 

820 

QA 

■ CINC) 


30 822 00 860 J = 1 , 2 

3 1 60 TO ( 83'1 |82H J I J 

32 . 82*1 NC o NB 

33 NO p 1^ 1 

3*1 IF IK) 824 , 826 1 032 

21 826 IF ^^C - 12> 830 , B28i 830 

36 828 QA « SV(NCil) - (SV(NCil) - QF ( l> ) • DT / I C ( 3 B ) SV(l8,lM 

27 GO TO 83*1 

38 830 QA ■ SVINCill {SVlNcTl) - C<^tn*DT/(C^3B) - SV(18,1J) 

1? 60 TO 83*1 

80 83 2 QA * SVINCiND) - (SV(nc,K) - SVlNCiNDH^DT / dITh ) ^ 

81 83** IJ r 2*1 + U - 2 

12 AB ■ VOL( I J ) 

13 AA » DT*«QA ♦ SV [NCiND ) J/2.0 

■*H 00 838' K s NO, *41 ^ 

HS IF (AA - AHI 836 , 836 » 8*10 

■*rr 836 AA m aA' * C I 36 1 • I SV 1 NC'i'K") =*-~S’V'(N’C , F + 1 D 

*17 838 continue 

"18 WRITE T5T8T01 I 

*»9 8*10 OA B AA- AB 

"STi K » 1C - ] 

51 IF (K) 882., 882, 8H4 

T2 gTl'2 OV « 5V ( NC 1 1) - QA — ^ 

53 IF (DV) 8B0, B*!*!! 8S0 

’5'M 88*1 DT ■ OA/(JA 

55 GO TO 860 

OV ■ SVlNCiK+11 ^ SV(NC,K) 








PRT*S RC9 







B6-f,a3*l32*TPFS.RC‘? 




-0!^ 


2 

3 


M 

5 

6 
7 

S 

9 

io“ 

t 1 


subroutine HC9 

niMENSION C(‘<0)i XNiR0,2», 5V(iai50), VlRAMlfl), RKIIRiR), 

J S C.I1H,5)» OCtlH)« A(6>» 0(15). F(2nl, VnL(lP). RP T 

2 QFUli TAU (5 J , CC ( 3 ) , CHBI3), CHl**), CPh(3) 

3 i:)<il M ) 


VlRANUfl), RK<lR,*n, 

■RI!iL(.2-L>- 


3 i:)<ii M ) 

COMMON/Z/ Cl XN, 5V, VTRAN* RKi SCi DC» At D, p, VOL, RPT 

j TAUt CC* CHH, Ch, CPHi piil t VE, Vl, CPflt CP.T* CADK i 

2 IRK, LOC, ITERX, INDEX, I, J, N, N 


» fiC , OF , 
_X 1. JT., 


2 IRK, LOC, ITERX, INDEX, I, J, N, N 

_C SETS VALUES I N VTRAN AR RAY 

CA949 FORMAIllH 7HSUB Rc9) 

DO 960 I = 1,5 


1 2 
1 3 


TA = TAUII) - (C(35) - 

LOC s T A7 D 4 1 *4 ) 

S V { 1 a , 1 ) ) 



I *4 
15 

902 

IF (LOC - H9 ) 9QR , 

WRITE (6,990) I ,LOC 

90H* 902 



1 6 
17 

90*4 

LOC a H9 
XLOC » LOC 




18 


TB a XL0C«0 1 1*4 ) 




1 9 


DT a TA - TB 




20 


GO TO (910»92Q»93D,9R0 

,950 » , I 



2 1 

9 10 

DO 9 1*4 J a 1,3 





22 

23 


29 

25 


26 

27 

~2 B~ 
29 


C lung to BRAIN C02»02,N2 TIME DELAYED ARTERIAL CONCENTRATIONS 
VTRAN(J) = KCF3 I J ) 


9t‘» continue 
C LUNG TO b rain 
VTRAN « 15) : 

GO TO 9 60 

920 DO 929 J s 
c brain to lung 


1* T 1 NE d elay ed 
RC> 3ri3') 


ARTERIAL concentration. 


*4,6 

CO 2 ,02 , N 2 


time delayEp venous concentrations. 




VTRAN ( J ) s RCF3 ( J ) 


31 92R CONTINUE 

32 ■ C BRAIN To'LUNG C OM'b I NED ' C02 , 0 2 » N 2 T 1 M E DELAYED VENOUS CONCENTRATIONS, 

33 VTRANI16) » RCF3(16) 

__ GO TO 960 

35 930 DO 93H J s 7,9 

36 C TISSUE TO LUNG C02,02,N2 TIME DELATED VENOUS c'c N C E N T R A f I 0 N S . ” ~ 

37 VTRANU) = KCF3tj) 

Tl 93R CONTINUE “ ‘ ' 

39 C TISSUE TO LUNG COOBINEO C02»02,N2 IImE DELAYED VENOUS CONCENTRATIONS. 

TO VTftAN(TT) = RCE3<r7"l ' 

*41 60 TO 960 

T2 VRO DO 9m*4 J = 1 , 3 

*43 C lung TO TISSUE C02,02 ,N2 TIME DELATED ARTERIAL CONCENTRATIONS. 

_ VTRAnTj + 9) = RCESrol 

*45 9RR CONTINUE 

*46 GO TO 960 

*47 C LUNG TO CAROTID SITE TIME DELAYED ARTERIAL CONCENTRATION. 

TB 9Ta VTRAN 1 I 3 ) ■a "RTH'3TTT) 

*49 C lung to carotid site 02 TIME DELAYED A RTERIAL TEN S J £N_» 

50 VtRAN(1*4)=RCH3(1*4? 

51 960 CONTINUE 

^2 C N ABEL I ST/OONM/'n R/\N 

53 return 

5*4 990 format (5X2ZHSV ARRAY ExTeEDED ON CYCLE 12,I2H WITH LOC a IS) 

55 END 







OTMENSIDN CMG)* XM‘tOi2)i S V > ©yV.ra^( I 8 J t RK(18»H), 

1 •5C(1H,5J» DCd'O* Ati», Dif5Vf^20), VoL(lQ)i RKT(2)i 


QF(A)» TAUtSJ, CCO), CHB(3), CH{H), CPh(3)f 

0 t H ) 


COMHON/Z/ Cl XN, SV, VTRANi RK, SCi DCi A, D, Fi VOU, RMTi BCi 
1 TAUi CC. CHB, CH, CPHi 0 Q , V£, VI, CPB, CPTi'CACK, X. D 


2 IRK, LOC, ITERX, INDEX, I| Oi Ki N 

C6969 FORMATUH HNSD8 RClO) 


: computes empirical functions fop acrdiac output and brain blood 

; FLOW differential equations 


C F(8) > PA(02I . 

IF (FI8 ) - 1 08.01 inns. I 02 f 1 . 1020 

C (DELTA )tl(n2) , EQUATION 7.3 . 

1 008 DO { I ) = ( ( -• 1 . 0U33E"5*F ( a ) + 2 ♦ 9 2 £ - 3 1 * F ( 8 ) .. D.28B5)«F(ai * 9.6AS 1 

C (DELTA ) QB (02 ) • EQUATION 7.9 . 

DQl2) = f((7.A5b«5E-b#F(«) - 2.32*)E’'S>*F(8) + 2.4(132E-3>*F(8) 


I - 0 . 1 323 ) *F C 8 ) + 2 . 785 

IJ^ ( D 0 I 1 ) ) l.P 1 2 » I Q I 6 * 1 G 1 6 

10 12 00(1'! = n.o 

ID 1 A IF ( DG ( 2 ) ) 1D29, 102S. 10 28 ' 

1020 00(1) * 0 .0 
1029 0Q(2> = 0.0 


C F(7) => PA(f.02) . 

ID2B IF (F(7) - 60.0 ) )(332 , I032 i 1036 

C 

C IF PC02 GT 60 I3QI3) STAYS AT ITS VALUE AT 60 - - OtD ROUTINE SETS 

“C THE V^aT tTE”"^” D N 1 3 ) EQUAL TO 0 

1032 IF (Pt7) - 90.0) 2036, 1090i 1090 


2036 D(3(3)s:Q» 
G0T01099 


C ( OELT A ) (5 ( C02 ) » REPLACES EQUATION 7.6 , 

1036 0Q(3)=6.D 


C 

GO TO 1099 


C ( DELTA ) Q ( C02 ) i EQUATION 7.5 • 
1090 DQ(3) = 0.3*(F(7) - 90.0) 


1099 IF (F(7) - 38.0) ID9B, 1052i 1052 

C ( DELTA ) (SB (C02 ) » EQUATION 7.11 • 


1098 D0(9) a ( 8 .(j 163E-9*F ( 7 ) - 3 . 1 07 3E - 2 ) * F ( 7 ) + ?. 3232E“2 
return 


- 99.0) 1056, 10S6i 1060 

1056 OQ ( 9 1 - 0.0 


Return 

C (DELTA)QB(C02) , EQUATION 7.13 . 
lT3“6 D"P'Q'r9'r"^nrTr~2TT7'9’8T^7‘rFl-7 ) * 9.391 S'E -5 ) •F'("7T'";r 1 ,299 7F“^l'* F( 7")' 
1 + D .7607 )*F17) - 15.58 


C NAMEL I 5T/D(a7DQ iF 

RETURN 


RT >S RCl 1 









original 1$ 
— OF-p0OR-Qofttnxr 


iba.(503*<32*TPF»«RC11 


SUBROUTINE- FtCll 


DIPENSION Cl'lDIt XMH0i2)| SVU8|50I, V7RAMl8»i RKIlS**!), 

1 SC(iSi5J, DtlH)! AU)i D(lb), F(20), VOL(IU), RMTI2), 


Rt(R), QFIA), TAlJ(5), CCt3), ChB(3>, CH(SJ, CPHC3)t 
DQ { S) 


, F . VOL, RNT , BC , OF , 

TAU» CCi CHB, CK, CPU, DO, VE, VI, CPB, CPTt'CADK, X, DT, 
IRK, LOC, ITERX, INDEX, 'H 0, P, N 
CALCULATES DIFFERENTIAL EQUATIONS 

1 ) 


EQUATION 10.] t 

ocin “ (vioctan - ve*c<i> ♦ Dt9)«(c(in»vTRANiH) + qFin 


I ♦VTRANI7) - C ( 10 3 *LC ( 1) ) I /C I 22 ) 

EQUATION 10,2 • ' ' 


OCI23 » (V!*C(32) - VE«C(?) ♦ 0 I 9 ) • 1 C ( 1 1 ) * V T R A N ( 5 ) + QF(l) 
1 •VTRAN(03 - C { 10 1 OF ( 9 3 ) 3 /C < 22 ) 


equation id. 3 • 



XABaS'5 • (ROT (23-. 215)46. -CUDS 

IFl {RP7(23 .GT.«2Jb)«AND,(XAB,GT»C,3 3DC( 10)»0c( ia)+AAfi/,030 


R2 

93 

C equation 7.7 . 

DC(ll) =* (-C(113 0,75 + DQ<2) + 0« ( H ) 3 7 C ( I 9 ) 


. 99 
95 

c equation IMD • 

DC ( 1 2 ) » F < I 9 3 / ( C { 39 3 *D ( I I 3 3 


>T5 

97 

C EQU AT I ON 1 < 1 1 , 

OC I 1 3 3 *= F ( 15 3 / < C ( 39 3 dD I 12 3 3 


Tff 



99 

TK 

DC«193 = F ( I 6) / 1 C ( 39 3 *0 ( 1 3 3 3 

- 


return 






BA_603‘432«TPFS.RC12 

1 subroutine RC12 

2 oniENSION C(HO)» XNI**0*2», SVMa.iSQ)* VTRAN(18)i RKUHf'U, 

3 I SC(l«t,5)t DC( HOr A-Ul . D<1S1« FUO). VflUlQ). RNT<7>, . 

H 2 BCmj, QFIA), TAU(5) , CC(3) , CHB(3), CHiRI, CPH(3), 

5 3 OQ ( H > 

6 CONMON/2/ Ct XN» SV, VTRAN» RK, 5C» DCi A, 0, Ft^'vOL, RMT, BC« QF» 

7 1 T AUt CC> CH8, Ch, CPH, OQ, V£, Vl, CPB« CPti'CAOK, X. DT . 

8 2 IRKi UOCi ITERX* INDEX, I, U, N, N 

? CQNMQN/R/ XDS>XMHiCXT.ftORK>OUMl >0Uh2.DUM3.ftQPK2iRMTB,RNTB2,T 1H EOF 

10 1 »RMLIN»ITTT 

1 1 COMHON/PLTBUF/NBUF,X8UF< leil ,^-BU'F ( 18I ,8) »NA(8) ,K5T0PP,TNAX,KPLT 

12 0ATANAyt PC02 P02 WqRK' 2 tH*) VI DB FReQAV02DF'/ 

1 3 DATA IRUN/*HUN * / i 1 ST OP / * S T 0 P « / , M OR E / * MQ R E » / 

t*! DATA I8ACK/*BACK»/ 

IS 0 I PENS I ON ftKKTT T ( SQ , 3 1 


!6 C&9A9 ■■ FORNATUH 8HSU8 RC12I 

17 C OUTPUT -- PUNCHED' CARDS AND PRINTED 


18 

CXT«C(35)+XDS-1D, 




1 9 

IFtCXT.LE.O* )CXT»+Q. 




20 

C DEAD SPACE VOLUME 




21 

DSVOL=O.IR0+D*OO2«VE 




22 

c respiratory frequency* 




23 

FREQ* «n.+ (*726»VE)/0SVQL)*‘**5-I.)/.363 





2M C DEAD SPACE VENTILATION 

IS DE ADVT = 1. «■» .098*V E 


26 C C(3U = (DEADV7*cn )+VE*0UM1 (DEADVT+VE J 

27 C C(32>°(nEADVT*CI2I'*VE»DUH2I/(DEA0v T * V E ) 

28 C C(33)=(DEADVT*C(3)+VE«DUN3>/rOEADVT+VE> 

29 C MNUTE volume* 

30 TVNToDEADVT+I VE+V 1 )/2. 

3j c heart rat e. 

32 HRATE®H3,8*(RMT{2I+C(26n+5‘I,5 

33 C 

3A C 

3S IFlCXT .LT. TIMEOF) GO TO 203 

T6 C I 

37 C HERE IF NEED TO READ A NEW WORK LOAO CARD* 

Tg C BRANCH IF IN BATCH MODE* 

39 TFMTTY .EQ* 0) GO TO 500 

TO C 


C 

Z here if TTY MODE 

IF< ITTT ,ED» 1 » GO 10 550 

C HERE IF TTY MODE, AND 1ST TIME THIS ROUTINE CALLED* 
ITTY - \ 

HR I TE ( 6,505 J 

505 FORMAT I ’ 0 I NPUT WORK CARDS****/ 

1 • ttORK= WORK LOAD (.watts I *«» */ 

2 * MINS= TIME FOR WORK LOAD*,**/ 

3 » PRINT- TINE INCRIMENTIMINSJFUR PR I NTOUT • • * * / 

R * EXEC. . . « / 

MOREb input MORE BEFORE EXEC,..*/ 

RUN B EXEC.WITh ABOVE, ThE-N CAN INPUT AGAIN,*.*/ 

7 * STOPb EXEC.WITH ABOVE THEN STOP*,.*/ 

8 * BACK* ERASE PREVIOUS WORK RECORD...*) 

50R ITTyIN *0 



2Q- 


- 

57 



ITTYOT « ) 

58 

C 



5T9 


HfH 

IF ( ItTY IT *LT . 50 ) GO TO 5D6 

60 

c 

HERE IF suffer FOR WORK LOAD CARDS IS FULL* 

61 



WRITEl6f51 1 ) 

62 


51 1 

FORMAT ( ’DBUFFER FOR WORK LOAD RECORDS FULL*’/ 

6 3 



1 • ft ILL USE EXEC« RUNr’ l 

6<» 



LEXEC = IRUN 

65 



GO TO 551 

66 

c 



67 


506 

ITTYJN X itt.yin + 1 

68 



IF (KSKP.GT.l .AND. KPLT.GT.O) CALL PAGE3 

69 


509 

VHR ITE t 6 *507 ) 



507 

FORMATl* WORK MINS PRINT EXEC •* 

71 



1 *(F6.2,1XiF6«2,1X,F6.2iIX»A8).»,’) 

72 



READt5,B02*ERR»509) (ttRKTTY(lTT'r,IN»0)*J»!*3),LEXEC 

73 


TOT 

F0RHAr(F6.?*lX»F6.2»lXiF6.2,IXtA8) 

78 



ttRlTE(6,503) (6 KKTTYUTTYIN,j),u«i,3),LEXEC 

75 


S03 

■ F'oR'M A T < 31* • *F6 * 2 ) 1 A8 » 

76 



IFILEXEC *NE. IBACK) QO TO 518 

77 



ITT YIN X ITIYIN - 1 

78 



IF (I tty IN .LT* 1) ITT YIN ■ 1 

79 



60 TO 509 

80 

c 



8 1 


516 

IFILEXEC .F.G1. IRUN *0R. LEXEC .EG. ISTOP) GO TO S5l 

82 



IFILEXEC .Ft|. MORE) GO TO 501 

63 



WRITE! 6,510) 

88 


510 

FORMAT!’ EXEC.PARAMETER WRUNG. TRY AGAIN.*) 

65 



GO TO 509 

86 

c 



87 

c 

HERE IF 1ST TIKE THIS ROUTINE CALLED* 

88 

c 

SEE 

IF more work CARDS IN BUFFER(WrXTTY!50,3)) 

69 


550 

IFIITTYOT .LE» ITTYIN) GOTO 55l 

90 

c 

HERE IE EXAUSTED ftORX CARD BUFFER i Vi R XT T Y ( 50 i 3 ) ) . 

CM 

c 

IptLEXEC *B.Q. IRUN) GO TO 508 
FORCE END OF COMPUTER RUN WHEN LEXEc- ’STOP*. 


^3 C( 15) = 0» 

GO TO 1210 


^ C 

96 551 W0RK2 » ftRlf IT Y ( IT TYOT , 1 ) 

^7 DURAT = V<RK1 T Y ( J TT TOT I 2 j 

98 C(39) = WRK.l T Y ( 1 TT.YOT , 3 1 












POOR QUALITY 


IF(WORK»SF. t5Q») TCT*2»3/(2.*«QRR/20Q.) 

IF (!»-nRK.l.T .50) TCTaH «6 


TISSUE 02 metabolic RATE. 

RMTf2)=SS0 2W(WQRK) -<SS02'Mft0RK)“RMTB2)«EXP (-TCT^ICXt -TIMEON)) 


VTlME=t»l-l.l*EXP(-TCT«fCXT-TIHE0K)/l,92) 
term used tM VI THAT 15 A COMPONENT OF TRANSIENT RESPONSE RELATED 


TO WORK LOAD. 

RMLIN =«SS02A{iftORK)-tSS02'A(ttORK)-RHTB2)«(l.-VTlME) 


IFTVTIHF. ,GE.l.) RMLINeSS0 2V)tWrORK) 
TISSUE C02 MFTABOUC RATE. 


RMT ( I ) = .8fl*RMT 5 2 ) 

IF(TVNT.GT.37*1 RMT(1 )»(TVNT-»R0«77)*RWT(2)/aa.5 


IFtC(35) .LT.CfRO) ) (jOT02 


IF (KP 


WRITE (6» 333 ) RHT ( 1 ) » R MT ( 2 ) 
FORMAT! * 0 * . I X , 2SHCHANGE IN 
1 5X,6HHR02» jFIO.H./) 


tiO TO 


! continu e 

U » AM0D(C(35) » Q.5) 

IF tU . LT. l .DE.S .OR. U ,&T. 
1F(C(35).LT«C(SO))QOT01230 
C (HO >«C ( SO ) +C( 3V ) 


arterial N2 TENSION. 


,H99V ) 


GO TO 1210 


1210 PAN2 * 

TISSUE 02 TENSION. 

PT02 = C(S»/.D(3) 

TISSUE N2 TENSION. " ' 

PTN2 » C ( 9 ) /D ( H ) 


10 PAN2 * D( 1 )*C(3) 


CEREBROSPINAL FLUID Ph , EQUATION 6.2 » 
PHCSF = 9 . - RCF 1 !Ch ( N ) ) 


VENOUS BRAIN H+ CONCENTRATION , EQUATION H»7 « 
HVB = CA0N*F t H ) / ( CC ( 2 ) - F ( N) ) 


VENOUS BRAIN PH , EQUAIION H.A . 
PHVfl » 9 . - RCF 1 ( h VB ) 


VENOUS TISSUE H+ CONCENTRATION « EtxUATlON 5,7 . 
HVT * CADKsP I 6 ) / (CC ( 3 ) - F(6J) 


VENOUS TISSUE Ph * EQUATION 5.6 • 
PHVT = 9 , - RCF 1 (HVT ) 


RESPIRATORY QUOTIENT (ALVEOLAR), 

RQ » ( ( C ( 1 I ) *VTR AN 1 H ) + Q F t I) • V T R A N ( 7 ) ) / C ( I C ) - CC(1))/ 


I « I^r9 ) » ICII1)*VTRAN15 ) + QF ( I ) * VT RAN(e))/C(lO)) 

QF ( 5 ) * QF ( 6 ) - RQ 


nMa(F{9)*CtlO>«.F(13}*(C(lO)-C{M))-F()2)*C(ll))*lOnO, 
AV020F = AV02bM/’C ( ! 0 > 









196 

T9’9 

200 




YBUn N6UF J »V I 
YOUr ( NBUF ,6 ) »C t 1 1 ) 


YPUF ( NflUF ,7 J «FKEy 
YB0F(NBUF,8)*AV02DF 


IF (kflt.gt«c> call plot 
IF (KPLT,GT*0) CjO to 1230 


HERE IF TTY OUTPUT* 


^RlTEl6,?on) CXT,RQ,QFtS) 

700 format ( *0T I h E * >F 10**11 *HINS * i 3 X » » AlV RQ*«F 1 0 * R . 


1 3X,*RQ DlFF*,FH,*t/ 

2 7X ,* alveolar arterial BRAIN TISSUE*, 6X, 

3 *C5F" V DRAIN V TISSUE*) 

ARITF, <6.701) C(1),CC(1),C(M),C(7),CC(2),CC<3) 


701 FORMAT!* C02 ' i R< F 9 , 8 ) , 9X i 2 ( F 9 • *t ) J 

WRITE<6,702) C<2)fFi9).Ct5>.C<e7,F(l2)»F{l3» 


02 ’ »‘UF9,‘n ,9X.2(F9*‘f) J 
03) C<3>,FtlO),C(6),C(9),C(6).C 


7D3 FORMAT!* N2 ' » 9 < F 9 , 8 ) , 9 X i 2 ( F 9 » 9 ) ) 

V1RITE<6,709) 0C( 1 ) »DCCR) .0Cl7) ,6c( 12) 


705 FORMAT!* JVAT • . F9 * 9 , 9X , 3 ( f9 . R) ) 

’■''iRlTEt6,7n6) DC(3)»DC!6)»nC(9),UC!lH) 


706 FORMAT!* IVES ♦ . F 9 , 1 , 9 X , 3 I F 9 , R ) I 
V1R1TE<6,7Q7) F(/),F(7).CPB.CPT,C(12),CPB,CPT 

707 FORMAT! * TTTiyi » . 7 1'F 9 , M ) ) ' “ 

IAR1TE(6,706> Ftl),F!niF!l7)»PT02,C(13).F(17),PT02 


709 FORMATS* OEM * i F-9 , 9 , 9 X , 3 ! F 9 , M ) ) 

WfR ITE I 6 ,705 ) DC!2) .DC!5) »DC(fl) ,L)C! 13) 


705 FORMAT! 
'■1R ITE t 6 


706 FORMAT! 
lAR ITE ( 6 

707 FORMAT! 
lAR ITE ( 6 


>AR ITE ( 6 

format ! 

WR ITE ( 6 


FORMAT ! 

lARITE < 6 


ORM AT [ 
■A R 1 T E ! 6 


WRTTE <6 
7 1 3 FORMAT ( 


1 

5X 


AR ITE < 6 
7)9 FORMAT! 


1 * • ,5X 

ARITE! 6 


,709) PAN2,PAN2fF!l8),PTN2.Cin),Fn8),PTN2 

♦ PN2 ■^TT(T9';9'j > 


7 10) !CH ( I ) , l3 1 ,9 ) ihVB ,HVT 


(H^) *,9)r»6!F9.9)) 

711) ICPh!t)tI=l,3),PHCSF,PHV8»PhVT 


PH * i9X i6 I F9 .9 ) ) 
712) (CUB! I ) »I«1 ,3) 


7 13) !TAU! I ) » I s lT5l , V 1 

transport times'/ 


,*AT*,7X,*AC*,7X,*VI*/ 

6 ! F 9 , 9 ) ) 


719) VE,C{10),C!1I),OC!IO)*DC(I1) 
**l2X*'VE*.8X»*Q*.7X,*Fe*.9x, 'DERIVATIVES'/ 


5 ! F9 ,9 ) ) 

7lB) FREQ,TVNT*OEADVT,HHATE,AV02DFiDSVOL 


1 , • DSVOL*/ 6X,6!F9,91 ) 


return 
















228 (^RlTE (*it9!0) CXTi RQi QF<5) 

229 C 

230 '»RITE (6,18151 (C(I), 1 » 1,3), 1DC(I), I » 1,3), F(7), F(l), 

23 1 1 ; 

232 VIRITE (6,1820) CC ( 1) , F(9), fUQ), Fl7), FM), pAtM2, CH(1), 

233 1 CPH ( 1 ) , CKB ( 1 ) 

239 write (6,1925) (C(l), I » 9,6), (DC{I), I > 9,6), CPB, F(17)» 

2 35 .J F ( 1 9 ) , CH ( 2 ) , CPH (2) 

236 WRITE (6,1830) (CU), I ■ 7,9), (DC(1), I = 7,9), CP T , P7C2, 

_2JU J °TN2, CH(3), CPH( 3J 

238 WRITE (6,1835) (DC(D, I » I2il9), (C(I), I = I2»l9), CH(9), 

239 1 PHCSF 

2B0 WRITE (6,18H0) CC(2), F(l2), C{6), CP8, F(17), F(l8), HVB, 

2jU I PHVB, CHB12) 

2*42 WRITE (6,1895) CC(3), F(131, C(9J, CPT, PT02, PT>^^, HVT, 

293 1 P HVT, CH8(3 1 

289 “ ’AR ITE ■ ( 6 , 1 850 ) (T'auII), I = 1,5), VT, VeT CdO), C ( 1 1 ) , DC (10), 

295 t DC ( 1 1 ) 

296 WRITE (6,1855 ) F R EQ , T V N T , OE AD VT , HR A TE , A V 0 2 D F , D S V 0 L 

_ 297 1 230 return 

299'” ' 1 290 format (5H XxX A 5 X 7 F I 0 . 9 ) 

299 1292 format (8F10,9) 

250 ~ 1905 format (iHl) 

25 1 1810 format ( 1 HQ6 X9HT IMEF 10 . 9 ,79X6MALV R Q F 1 0 , 9 . 3 X7 H R Q D1FF«F8.9/ 

252 1 16X3HC028X2H028X2HN27X21HO E « I V A T I V E S9X9HPC026X 

25 3 2 3 H P 02 7X3HPN27 X 9H ( H> ) 7 X2 HPMSX9HH0O2) 

259 1815 format (3X8HALVE0La'r9F1D»9) 

25 5 1820 format (3X8HAR 7ERTAL3FlO«9,30X,5F1 0»9,Fa.9) 

256 1825 FORMAT ( 6 X 5h B R A I N 1 1 F 1 □ , 9 ) 

257 1830 format ( 5X6HT I SSUE 11 F 1 0 . 9 ) 

25’S TT35 format (8X3hCSF3QXBF 10,9) 

759_ 189 0 forma t (9X7hV BRA IN3F10,9,3QX,5F10t9,F8.9) 

260 1895 format (3X9HV T 1 S SU E 3 F I Q , 9 i 30 X , 5 F I 0 • 9 , F 8 , 9 ) 

261 1850 format ( 6x IShTRANSPORT TIMES --9 X 2 H A B 8 X 2 H V 8 8 x 2 H V T 8 X 2 H A I 8 x 2 h A C 2 X 

T5T' ■ i X'2 H VlTrrHWe*TTFQ 9T^FB‘7 'xT^ I V aTTv rf /'2 1 X , 1 0 FfoT 9', F 8‘r9 ) 

263 1855 FORM at ( 3X , 9HRESP F R E Q , F 8 . 9 * 2 X , 1 3 H M 1 N 0 T E VOLUME, F8, 9, 

j 2X,8H0 5 vent ,F8 ,9 ,2X , IOHHEAR r RATE,F8,9» 

265 2 2X ,7HAV02DF ,F8.9,2X ,5HDSV0L,P8.9 ) 

?T6' C ‘batch MODE“wORK CARO READ.,, ' 

267 C 

“ 26l C WILL U5F WORK CARD WITH T 1 M E ^ A S ‘l D I C A tVo‘n ' ' 

269 C OF END OF RUN BECAUSE 1)06 HAS PROBLEM 

?713 t WITH Ends On read, 

?71 500 REAO{5,300»END=»2) W0RK2,DURAT 

27T Jffo FWma'T ( fTT 2T3X"» F 6 i'2') ‘ 

273 C 

T7"9 IF t OURAT 70T, QO ^ TO 606 

275 C HERE IF READ INDICATION OF END OF RUN IN BATCH MODE, 

Tfl, C ( 15) = 0. 

277 (30 TO 1210 

278' C ■ 

279 606 IF (KPLT,GT»0) CiO TO 607 

280 write ('6,305) WORK 2 , D UR A T , Cx T 

281 ■ 305 FORMAK ’0* *93 ( 4* * ) / 

282 I * WORK LOAD ChO . ( • , F 6 , 2 , • W AT T 5 FOR*, 

283 2 F8«2, * MINS) AT ’ ,F9,9 , *M1 N5» 1 _ _ 

289 607 T IHF0P = '7UraT + CXt' 




-0RfelNAfc-Fji^(j£-;?a 
POOR OUAU ^ 


time constakts for work loads 


TISSUE 


SYSTEM responses 


metabolic rate* 

IFItt0RK2.GE»W0RKJRMTB2«RMTt2) 


CREASING ttOR 
IF {ftORK2 .LT •WORK ) RMTMaRMT I 2 ) 


0RK2,l.^'>W0RK)RfnBaSS02»lW0RK2) 
IFI(WORK2.LT.WORK).ANO,(WORK.GE*5Q.U TCT = 2»3/(2.*W'ORK/20Q.) 


IFMW0RK2.lt. WORK). AND. (WORK, LT» 5Q#))TCT=S*6 
1F(WORK2,GE.WORK1 GOTOl 


01 W0RKbWoRK2 

MARKERgQ 

TISSUE 02 METABOLIC RATE. 

RMT(2)=RMTB-{R«TB-RMTM)*EXF«-TCT*(CXT-TIME0M*.50) 


VT1ME=1.I"1»1*E)'P(-'TCT*(CXT-TIME0N)/3.8R) 
term used in VI THAT IS A COMPONENT OF TRANSIENT RESPONSE RELATED 
TO WORK LOAD. 

RMLIN «RMTn-(RMTB-RMTM)*(l,-VTlME) 


F ( VT 1 ME »6E . 1 . ) RMLINcRMTB 
TISSUE C02 metabolic RATE. 


RMT( 1 )b.88*RMT(2) 

IFtTVNT.GT.37.) RMT!l)»(TVNT+R0.77)*RMT(2)/efl.5 


IF ( C < 35 ) ,LT»C t MO ) ) G0T02 
IF (KPLT .GT «0 ) GO TO 2 


WRITE (6»333) RMT { I ; »RHT ( 2 ) 
G0T02 







iBA.,G03*432<iTPFS*Rc13 



SUBROUTINE RCt3 


dimension CISO)! XN«H0*2J, Svn8,5D>» VTRAN(I8), R|f(lH»‘n, 

I SC(l«i»B)» DC(lR)i OU5)» F(2p)t VoHlQ)t RNT(2) 


RCIR)* QF(6), TAU(5)i CC(3)i CHB(3), CHI')!* CPh<3), 

0 « ( 4 ) 


SOLVES M differential EQUATIONS By FOURTH-ORDER RUN&E-KUTTA AND 
AD AMS-MOUTI.ON M RE D 1 C TO R -C 0 R R ECT 0 R METHODS 


NAMELI5T/DHQ/C»DC,SC 


COHMON/Z/ Ct aN» 5V. VTRAN» RK, SC« DC* A, Di F» VCL, RMT* BC. QF* 

1 TAUi CC» CH6» Ch, CPH, DQ < VE, VI, CPR , CPT*'caDK. X. DT. 

2 ' IRK. LOC, ITERX, INDEX, J, 0, M, N 

FORMAT i IH 8HSUB RC 13 ) 


IF ( IRK - *»J 13( 


1304 DO 13S2 INDEX » 1,4 
DO 1308 I » 1,M 


RK n ♦ 1 NDEX ) s DC ( I) 
1308 CONTINUE 


GO TO < 1 3 1 2 , 1 320 , 1 , 
1 3 12 DO 1 3 I A I 3 1 ,M 

SC n , I RK+ 1 ) 3 ecu 
SC II , IRK » = DC I M 


1316 CONTINUE 

TI » COS) 


1328, 134o>, index 


24 

25 

1320 

CI35) 8 TI + CO6)/2,0 
DO 1324 I 8 1,M 

26 

27 

1324 

CII) « SC(IiIRK+l) + C 06 ) *RK ( 1 , INDEX ) /2 ,0 
CONTINUE 

29 

1 328 

GO TO 1336 

COS) 8 TI ■> C06) 

30 
3 1 


DO 1 332 I a 1 ,M 

cm * SC(I,IRK+l) + C06)*RKI I »INDEX) 

32 

33 

1332 
1 336 

CONTINUE 
CALL RC14 

34 


60 TO 1352 

35 

1340 

DO 1384 I a 1 ,M 

36 


cm « SCtI,lRK*l> + C06)*IRK( 1 ,1 ) ♦ 2»0*RK(I,2) + 2.0*RK(I,3) 

37 


1 + RK ( I ,4 ) ) /6 »D 

3 8 

1344 

continue 

3 9 


IRK a IRK + 1 

40 

1352 

CONTINUE 

41 


RETURN 


= I * 
SC t 1 ,B ) a C < I ) 


SC I I ,4 ) 8 OC ( 1 ) 

cm = sell, 5 ) + C06)* (SS.O*SCl I ,*♦) - 59,Q«SCn,3) + 37,Q*SCU,2 


1 - 9,0*SCI I , n 1/24.0 

1360 CONTINUE 


C ( 35 ) = COS) + C ( 
NC35*C ( 35 ) /C 06) + .1 


C(35)»C 06 ) •NC3b 
1364 CALL RC14 


00 1368 I 8 1 ,M 
SC n , 1 ) « C ( I) 


cm » sen, 5) ♦ Ct36)Fl9.Q*0Cm + 19,0»SC( I , 4 ) - S.0«SC(I,3) 
1 + sc ( I ,2 ) ) /24 .0 

068 continue 






S7 

DO 

1 37 2 I 

» I ,M 



58 

IF 

( AB5 

1 C (1 ) ^ SC t I 1 1 ) } - 

1 ,0E-3) 

1372» I372f 136^ 


1372 continue 

00 1376 I a i,M 
075 1 356 J » 1,3 
SCI 1 lO) » SC (I » 0+1 I 
3 1376 continue 

^1 RI-TURN 

5 END 




COHMON/Z/ C* XNi SV, V7RAN» RK, SCf DC» A, D, F| VOL, RMT, BC, QF, 

! TAUt CC» CHS. CH, CPH. D lS| i VE» Vl. CP,B„t CPT » CAD K . X . DT ♦ 

2 IRK. LOC, n£RX. INDEX, I, O, H. N 

c CAULS ot her sobkoutines in a block _ 

10 C&969 FQNMATtlH aHSOB RCIR) 

U ALL R C 3 






fv CO o jo m 





nrs^6"G a'oa^TPFs.RciB 


j*AeE_i^ 

QUALTO 


I 


SUBROUTINE RCIB 

"imiTR'STDir rnfo ), svufiisoj. vTra vri'e T,~^ 

1 SC(1H,5»» DCIM)* At&N D(15>» F(2a)i VcLllDl* RNT(2). 

2 " "Birni") . QF'( 6) ,”'T4 'ot 5 ) ; cco)', c BTyyr c h ( m j 7 cPHTaTT 

3 015(H) 

COr^MON/Z/ Cl ^Ni SV , VTRAV t RK , SC » DC « A , Di F, VOL, RKT» BC. <5F. 

,_L_ T A{j I _cc_i CHB » _.C H j CPH t 0(5 i V V. I •_ C P_0 » C P T » C K j _i! i r. T , 

LOCj ITERXi INDEXi Ii J> «» N ' 


IRK 


C6969 


FORMAT ( 1 H 8HSU8 RC 1 5 ) 


C NAMEL 1 ST/SCh/S V 

C SHIRTS VALUES IN SV ARRAT 


IH 
t 5 
TS" 
1 7 

*rg- 

1 9 


DO 1 S3D I = 1 , IM 
DO 1520 J =■ ItHV 


JM = 5 1 - 0 

JMM a JM - 1 

S V t I Tj ^ ’) ^“S'V "( T» jmmt 
1520 COMTINUF 
1530 CONTIN U E 


RETURN 






I 5V M , 1 ) » CC ( 1 ) 

15 C arterial 02 C0^CEN^RAT10 ^» 

16 SV<2,1)=F19) 

12 C BRAIN VENOUS C02 C ONCF.N T R A T. I QN , 


1 

' c 

sv ( 4 s n 
arterial N2 

= CC (2 » 

CONCENTRAT ION. 



! 

sv n , 1 ) 

= n I Q ) 



C 

BRAIN VENOUS 02 CONCENTRATION 



i 

SV I 5 , 1 ) 

= F( 12) 



c 

BRAIN VENOUS N2 CONCENTRATION* 



1 

sv 16 , n 

» C ( 6 ) 



1 c 

TISSUE VENOUS CQ2 CONCENTRATION, 




SV ( 7 , 1 ) 

» - CC ( 3 ) 



c 

TISSUE VENOUS 02 CONCENTRATION* 




SV < fl , 1 ) 

= Fill) 



c 

TISSUE VENOUS N2 CONCENTRATION* 




SV ( 9 , I ) 

o C ( 9 ) 



c 

CARDIAC OUTPUT* 




SV ( 10 . n 

= c ( 1 n ) 



c 

cerebral BLOOn FLOW. 




SV < 1 } , 1) 

a C ( 1 1 ) 



c 

TISSUE BLOOD FLOW. 




SV( 12, IJ 

a QF ( U 



c 

arterial h+ 

concentration. 




SV( 13. 1 ) 

= Ch ( 1 ) 



' c 

arterial 02 

■tension. 



t 

SV ( M , 1 ) 

a F { 1 ) 



c 

initial t ine. 



t 

SV (15 , 1) 

a 0,0 



1 c 

total gas concentrations at brain 

EXIT. 



SV( 16. 1 > 

a SV(H.i) * SV(5,n ♦ 

5V ( 6 , 

1 ) 

> c 

TOTAL GAS CONCENTRATIONS AT TISSUE 

Exn 

• 


SV( 17. 1 > 

a SV(7,1) + SV(8,I) * 

SV(9, 

1 ) 

c 

simulated time. 



1 

SV ( 1 8 . 1 1 

a C { 3B ) 



\ 

return 















5 


SUBROUTINE «C17 


1 

“2 

3 

5 

“S' 

7 

9 
TO 
t ! 
TT 
I 3 
19 
1 5 

T& 

17 
1 8 
1 9 

I 


I 

r 


I 


35 

37 
TB" 
39 
“911 
9 1 


13 


ITTMEWrON r.l90>i XM9072), bVllB»50<, VTRAN(i 8 )t RKU 9 i 9 ), 

1 5C(19,5)» DC119)» A( 6 )» 0(15), Et20)» VQL(lO)t RMT(2», 

2 "^ BT'{ 9 ) , U'nTDT'^TA u ( 5 ) , c Cl 3 ) , C H B t 3 ) , C H ( 9 J , C P H ( 3)“, 

3 0 Q ( 9 » 

CTTFFTfNTTT/ C, XN, ST;~TTirA N 1 Rk. SC» DC. A, D, F, VOL. .RHT, BC, QF , 

1 TAUi CC, CHS, CH, CPH, DQ, VE, VI, CPB, CPT» CADK, X, DT, 

2 IRK, UOC, TTer jT7 Inoex,i,.j,h«n 

OOMMON/R/ XU5,XMh,CXT,A0RKiDijMl»DUM2,DUN3,W0RK2,RNTU,RPTB2,TIME0F 
j ,RmL fN 

C NAMELIST/BAD/CH19) ,CADK,0(ll ) ,C112> ,BC(9J ,C(37J ,C{30) ,VTRAN(19) , 

C !TFRM,VI,C{2Cl),C(16),VTRAN(15),Cl2P,VTRANn3),C(37)»D{9),Ctin, 

C 2VTRAN{U),qFU),VTRAN(17),(.({0),F(ll), 

C6969 FDRMATIIH 8HSUB RCI7) 

C CALCULATES VENTILATION 

C CFS H+ CQNCENTRA TION , EQUa’tION 6ql • 

CH<9) = CA0K«D 1 1 1 ) «C { 1 2 ) /BC ( 9 ) 

IF (C(37) ,&T» l.OE-5) <30 lO 1708 

1709 VI » C(38) 

50 TO 1730 

17Q8 TERM = 0,0 

C DECISION '0N“ ARTERIAL 02 TENSION AT CAROTID BODIES'SITE* 

IF IVTRAN(19) - 109 »D) 1710, 1720 , 1720 

1710 term = (23. 6E-9)*( < 109.0 - V T R A N C 1 9 ) ) * « 9 , 9 ) 

C CONTROLLER EQUATION AS A FUNCTION OF HUMORAL TERMS* 

l'7 2 0 VI s rr 2 0 ) ■• 1 C ( I 6 ) *VTRAN{15» + (1,0 - C(16))*CH(9)) 

1 + C(21 )*VIRAN( 13) + TERM C(37) 

"C INCLUSION OF NEURAL COMPONENT AS A FUNCTION OF ftORK L(3AD* 

SVNT 2 = SS VENT ( SS02Vi t VtORK ) ) -VI 

IF( (SVNT2.GT.O,».ANO.(SVNT2.LF. ,15.) > VI»V I+SVNT2 
IF(SVNT2,6T» 15, ) V1=VI+15» 

J, - - 

C DESCRIPTION OF TRANSIENT VENTILATION RESPONSE* 

SVNT aSSTTN'TTRMLTN. ) ; 

IFISVNT.GT.O.S) Vl«VI+0*7S*SVNT 


C EXPIRED VENTILATION RATE, EQUATION ll.l , 

I 7 30 VE'"*"^! ■ + 0 I 9 ) * rc ( I I ) *VTR AN l I 6 ) + QF < 1 ) <* VTRANTrT) - C(10*F(nT) 
IF (VI *LT, 0*0 .OR, VE ,LT. O.Q) (50 TO 1790 

RnnwTv 

1790 VI a 0.0 
Q--TJ- 

return 



















)B4-S03*>32*TPFS«RCI9 

1 SUaROUTINE RCl^ <CPA> CVHBA, CVCt BHrAt fc) ^ 

dimension CIHD)* XN(Mo»2). SV(l8fS0)t VTRANU8)* RKtlHt*»)i 

SC(lH,5)t Ocn*ni AU>> Dtis)t F(2 o). VOUlQl. RMT(2). 

SC(R), QFtA), TAO(SI, CCt3), CHBI3), CH(R), CPH(.3)t 

J DQ I R > . 

C0MMf)N/2/ Cl XN, SV, VTRANi RK , SC* DC* A, D, Fi VOL, RMT| BC| QF, 

J TAUt CC< CHB* CH , CPH, DQ, VE, VI. CPR . CPT»~CADK. X, DT, 

2 IRKf LOCi ITERXi INOEXi.Ii Oi M, N 

NAWELIST/DM2/CPA,CVhBAiCVC»8HCAtFC 

10 C6969 FOHMAHIH 8HSU8 RC19) 

JJ C iterates FOR venous brain and V&N005 TISSUE CQ2 CONCENTRATION 

12 C term used in equation *f,2 • 

13 1910 X »» (CVC - FC)/<0»01«CPA) 

19 c logarithm subroutine. 

}S X » RCFl t X ) ^ ^ 

16 C EQUATION 9. 2 . 

17 X » BHCA + Q,37S»CCtl7) - CVHBAl - Die)«(X » 0.19) * FC 

18 CALL RC6 {CVC> 

19 CVC B CVC * 2.D«(X « CVCl/3 .0 _ 

20 IF UTERXI 1920. IrTo, 1920 

2J 1920 CONTINUE 

22 return 

23 EM 






ORIGINAL PAGE I& 

JBA-6b3 9 3Z*TPFS .Rc20 OF POOR-QUAUW 

1 SUBROUTINE RC20 

DIMENSION t«Hbf » XN«*»0i2)| SVIt8i50»f VTRANH8)| RK»l«li‘t), 

SC(M,5)i DC(1**)» Al6)» DU5)» F(20)» VoLClOJ* RKT(2)» 
BCI*ni QF46)i TAUIBI, CCt3I, CH8I3)* ChlSit CPHI3)* 

I OS(HJ 

i COMMON^Z/ C» XNt SVf VTRAN* SC* DC* A| D| Ft VOi-i RMT« BC • QF> 

[ 1 tAU* CC* CHS* CH, CPH* DQi VEt VI, CPB* CPT * ' C ADK i X , DT, 

} 2 IRK* LQCi ITERX* INDEX* I* M> N 

? C NANELIST/NMF/F 

C6969 FORMA7UH 8HSUB RC20) 

C SETS TIME DEPENDENT EXPRESSIONS 

12 C arterial oxygen CONCENTRATION INCLUDING EFFECTS OF HEMOGLOBIN* 

1 3 F(9) » 0t6>»C(2) ♦ CHBM) ^ 

IR c arterial nitrogen concentration, 

IS FMO) « D(7)«c«3) 

T5 c nmn: Artemi A L gas concentration at lung exit* 

17 F ( 1 n ° CC (I ) * F ( 9 > •» F ( 1 01 

TS C VENOUS brain oxygen CONCENTRATION INCLUDING EFFECTS OF HENOGLOBIN* 

19 FU21 = CIS) ♦ CHBI2) 

TO C VENOUS TISSUE OxYGEN CONCENTRATION INCLUDING EFFECTS OF “ hEHQGLC B I N * 

21 FII3)=C(S)* ChB ( 3 1 

E oxygen tension IN BRAIN. 

23 F{ 17) » C(5)/D(3) 

t NITROGEN TENSION IN BRAIN* ^ " 

25 F I 18) « C I 6 ) /D M ) 

C PRODUCT OF diffusion COEFFS.AND D IFFERENT 1 Al“TRaIN - CSF GAS TENSIONS 
FIH) «= C(27)«( CPB - C I I 2) ) 

F(15)»C(2e)*lFtI7)-iC(l3)) 

FM6) =■ C(29)*(FI18) - CI1*M) 

C 


RETURN 

END 
















iB6-G03*)32*TPF#oRc2 1 
1 5U 



C49&9 


SUBROUTINE RC21 ICHBA, FA, FD, CCAi CHA, CPHA) 


DIMENSION C(‘iOJ» XMMO,2>, SV(ie,50»t VTRAN(18), RKUHt*!), 

J SCU**»6)t DCIlRIi A < ♦ D US > , Fl20)f VoLUOJ* RMTI2) 


BC(H), QF(6), TAUIEJ, CCI3), CHB(3), CH(HJ, CPHI3), 
OR ( H ) 


COMMoN/2/ C» XNi S.V, VTRANi RK, SC» DCi A, D, F» VOU, RMT, BCi i 
TAU, CC, CHB. ChjCPH, 


2 IRKi LOC. ITERX, index* I* J> Hi N 


69 FOhHATUH 8HSUB Rc2l) 


NAMELlST/PR/CHBAiFA,FO,CCAtCHA,CPHA 
COMPUTES h+ lONi Ph, and OXYHEMOGLOBIN 


arterial h-»> concentration. 


arterial 

CPHA > 


PH. 

9.0 


- RCF 1 ( CHA) 


U 
) 7 
1 8 
1 9 


: DEVELOPMENT OF EXPRESSION USED IN CALCULATION OF ARTERIAL 

C OX YHFHOGLOB I N SATURATION. 

X » RCF2ICPHA) 

X a • fa 



SPRTiS RCF) 












RAGE 

OF POOR QUALITY 


3B6-603*»32oTPFS«RCFI 


1 


FUNCT ION RCF 1 f I 


2 

3 

~>r 


OlMFNSrON 


J 


CtROJi XNtRQ,2)f 
SC ( 1 *♦ |5 } . DC ( I'D t 


B C t H I. » 
0 HI R 1 


SV(iHi50)» V.TRAM18)* RKllR.R), 

A(6)i 0(15J, F(20l« VoUlCIi RPT(2), 
filFl6), TAO{5)i CCI3), CheOli CH(HJ, CPhO>» 


COMMON /Z/ C» XNi .SV, 

TAM* rr . CHB I 


VTRAN » 
CH, CPHi 


D, P, VOL, 


6 

7 

T" 


T AU » 
“~rRl <7 


■ 9 

Tier 


1 1 


cc , 

L'^'i ITERXi 
BASE 10 


RK, SC» DCi 


I 

DQ , 


V£i 


jtTN, uu*., iirnAi iNutAi 

logarithm to base 10 

RCF 1 a 0*H3‘T29‘<Ma * ALOMt^} 

return 


VI , 

INUEx; rr"j| M* N 


- RMT , ec , 
CPS, CPTt CADK, X, DT, 


QF 





OF POOR 


)BA-ri03‘<32*TPFS ,rcF2 ' 

\ ruNCT 1 ON RCF2(Z) 

DIMFiNSION CMni. XNtH0.2)i SVUB'jSO'), VTRAN(18)t R K C 1 » ‘t’) , 

SC(lHt5)> PCtiH). A(6)t DUS)» P(2rHi VnUlDli 

BCt**!* QF(6>f TAU(S)t C C < 3 ) , C HB ( 3 ) , Ch(*U, CPH(3)t 

! DQ ( 

COhMON/Z/ C» XM, SV, VTRANf RK, SC» DCi A, D, F, VOL, RPT, B,C Ci'F t 

1 Al, I C C I CHBt Ch, CPhi Dtii VEt VI, CPB, CPT» FADK , X . DT . 

2 IRKi LOCi I7F.RXf INOEXf 1» o$ Ni N 

C OXTHEN ri^L 0^1 N - PH EPP 1 RICA L FUNCTION ; 

c equation i.”'* . ' " ' 

RCF2 = tlta»00&68J5*Z> *" D»1UC90)*2 + - 0«R5R 

12 RETURN 

13 END 


?PRT*S RCF3 


1 n 
1 1 


1 



FUNCTION R C F 3 { K K > 




N5I0N Cl-IDJ* XNt‘»Oi2J, SV(|Bi50»i VTRAN(18)» RKtlR»H), 

SCCIH.SJ, OCIlHJi A<6>, OtlS), F(20). VqUIO). RFT{2), 


? TAU(5), CC(3), CHBUJ. Chi'll, CPH(3), 

1 Ofil'O 


I FiiDJ, VQLVIOJ, Rf'TI 


COMmoN/Z/ C» XN, SV, VTRAN* RK, SCi DC» a, 0, F, VOL, RWT, BC 
1 Taui CC • CH8, CH, CPHi DU, VE, VI, CPB» CPT« CAOK, X, 


2 ,IRK» UnC, ITERX, index, 1, U. P, N 

IJ « A « •! A » I ^ ^ I A 


VTRAN function 


lABLES WITH TI 
RCF3 = SV{KN,L 



'PRTtS SS02'-V 




Toy 


5B6^603M32*TPFS«SS02ft 

1 FUNCTI0NSS02’A(XJ 

: calculation of steady-state OAYQEN REUU IREHENTS for VARrOUS levels 

C OF WORK LOAD (X«ftATT5)» 

IF C Y ,GT .250 9 )G0T01 
IF( X .LT.75. > G0T02 
SS02fta-»U72+tX/70» ) 

Return ^ 

l SS02i\«3.5 

RETtiRN 

ID 2 SS021A»I X/75« ) ♦.2 15* ( 75.-X ) /75 . 

lie 

12 RETURN 

1 3 END 









U 


T 

FUNCTION SSVENT«X> 

C calculation of STEADy-STATE VENTILATION RATE AS A FUNCTION OF TISSUE 
C OXYGEN METABOLIC RATE. 

IF < X ,Lt . .2 Ih ) bbVTNrcSiSTR 

!F(«X»GT. .215). AND.IX.lt. 2. ))SSVENT“25.*X 


■5 



5B4«GQ3H32*TPFS.PtOT 

1 SUBROUTINE PLOT 


2 

3 


COMMON/PLTBUF/IPLTPT(TLOTBFUaiI»pLOT0FI18lt8)iALPHA(8)»KSTCPP 
(, »TMAX,KPLT 


R 

5 


dimension HEADI6) 

DATA HE AD/ * 6R0D 1 NS RESPIRATORY CONTROL MODEL 


6 

7 


DIMENSION TSTEPA I R ) »TSTALP I *♦ ) 
DATA TSTEPX/1 « >1 » ,60o, lRRn«/ 


8 

9 


DATA TSTALP/*KINS**»MINS»,»hOUR»,*OATS*/ 
dimension II DAS( 8 , 2 I iPPARS < 8 ,3) 


I 0 
I 1 


dimension XNO(B) 

DATA I5A/«S VtIN/«N 


*/. IT/'T «/tIPLL/*PLOTt/ 


12 

13 


DATA IE IRST iKSTOPP ,N0P/0 .0 »D/ 

IF (IFlRST.EQtO) CAUL IMTT(300> 


18 

15 


16 

17 


K»8 

IFIRST=1 


1EXECI»1PLI- 

IS B IPLTPT 


• I 


1 8 
1 9 


20 

21 


IFIIPLTPT , 
iFt l EXECI . 
1 GO TO 28 
IFtlEXECI . 


.LE. 


1 ) 60 
IPLL 


TO 1 
• AND. 


XST OP.P 1.IQ •. 

£G» IPLL ,AND» KSTOPP 


0 ) 


6E« 1 ) 


22 

23 

28 

25 

26 
27 


1 6o 
IS 


TO 

I 


I 20 


1 

II ftRITEl6.2) 

2 FQRMAT(9X|* graphic uk;irui \ i « w »J3 ; jL i._» n 
6 I ♦ ST OPX-, ( A2 , 3F5 .0 J « . • ) 

8 READ(5«3iERRb 8 > I » TSPL I t *STARTX tSTOP X 

3 FORMAT! A2,3F6. 01 

IFtI »EQ« IN) RETURN . 


28 

29 


OUTPUT ( Y .N » S ) xl 1 M_E U' T E R VA L S } S T ^ T X » • 

« . • 1 


3Q 
3 1 


32 
3 3 


38 

35 


IF < I • EG , 15A ) 

IF CSTOPX»LT.t 


GO TO 10 
i> STOPX»TMAX 


IF (STARTX.lt* l» 1 ST ART XaTLOTBF ( 1 1 

RUNSTPsSTQPX.STARTX 


TTSPT = TSPLIT 

.IFtI *NE . IT) GO TO U 

iFtTSPLIT «LT» 3. ,0R. TSPLIT .GT. S.) GO TO M 


HERE TO aUlLD PLOT PARAMETERS. 


36 

37 


38 


NOP a 0 


39 


00 20 I a 1 1 K 


80 

17 

VlRITEt6,5) ALPHA! 1 ) ^ 


81 

5 

FORHATt* •,A6,13X,»T SCALE ( A 8 » 8 X , F 8 . 0 , 2F 6 , 0 ) * / 


82 


] » plot ( Y |N *5 1 LOC HIGH * , . ’ ) 


83 

7 

READ ( 5 , ^ »ERR“7 1 I 1 »XPI »XP2iXP3 


88 

^ 

FORMA’! ( A8 ,8X ,F8 .0 ,2F 6.0 1 


85 


ipt 1 1 .EQ. IT ) GO TO 18 


R'S 


I F t I 1 .EG. 1 S A 1 “G 0"T0"T'9” 


87 


if! I 1 .NE. IN) GO TO 17 


TB 


FPARS < I f 1 1 “ 0 . 


89 


GO TO 20 


50 

18 

IFtXPI .LE. 0. ) GO TO 17 


51 


IF(IFlXtXPl) .GT. 8) GO TO 17 


^ 


IF!XP3 .GE. XP2) GO TO 17 


53 


PPARS t I , 1 ) = XPl 


M8 


PPARS ( 1 , 2Y o'“5rP*Z 


55 


PPARS ! 1 ,3) a XP-3 



56 


19 IFt IFiXtPPARSt I til ) .GT 


NOP ) NOP a PP ARS ( I » 1 ) 















5 7 1 F J PP ARS (III) nrri TT) GO TO 17 

58 20 CONTINUE 


■ 5)9 
60 

C“ 

BUl'L'O VlihGUF. page iRAPh,; 
10 CALL NEHiPAQ 



1 

^2 



CALI. 

CALL 

MOVAbS ( 300.77 5 ) 
DNPBUF 



63 



ViR 1 TE 

16,21) HEAD 



69 


2 1 

format (25X .6A6 ) 



65 



CALL 

M0VAB5(3,1D0) 



66 



CALL 

DR^^Aes 1 lonci , 1 00 ) 



67 



CALL 

DR'^ABS { 1000,750 ) 



68 



CALL 

DRftABS ( 3 ,750 ) 



69 



CALL 

ORttABS (3,100) 



70 



call 

MOVABS ( 300,750 ) 



7 I 



CALL 

0RWA8S ( 300 , 100 ) 



72 



CALL 

DMPBUF 



73 



M = 

mci JM 



79 



12 = 

300 



75 



1 3 * 

700 / (11 ♦ n 



76 



DO 25 

I “ 1,11 



77 



12 = 

12 + 13 



78 



CALL 

MOVABSl 12,100) 



79 



TTLL 

"DR-AABS 1 12 ,750 ) 



80 


25 

CONTINUE 



81 



CALL 

D M P 0 U F 



82 

c 

DRA« DIVISIONS FOR DIFFERENT GRAPHS, 



H3 



137 e 

650 / NOP 



89 



I 9 Y = 

I3Y / 2 



8S 



I 1 = 

NOP - 1 



86 



12 = 

100 



8 7 



I F ( I 1 

♦LT . 1 ) GO TO 280 



88 


- 

DO 26 

I = 1,11 



89 



19 = 

12 + I9Y 



90 



CALL 

NOVABS ( 300 ,19) 



9 1 



CALL 




92 



CALL 

MOVABS ( 995 ,19) 



93 



TTtnr“ 




99 



12 = 

12 + 137 



95 



CALL 

MOVABS C3 , I 2 ) 



96 



CALL 

DRWABS (1000,12) 



97 



Call 

WFHUF 



98 


26 

CONTI NUE 



9"9 


^ a 0 


1 1 + T9-T 



1 00 



CALL 

MOVABSOOO, 19) 



rm 



"cTTrcr' 

ONVCABSISOS, 19) 



102 



CALL 

MOVABS (995 ,19 ) 



ftnj 



CALL 




109 



CALL 

DMPBUF 



prrs 

t 


1 A UPH A 



106 



I ! = 

750 



1 U7 



IDUP 

» 1 



108 



CALL 

VftINDOlO., 1023*, 0., 780*1 



nr? 



CALL 

SWINDOIQ, 1023,0,780) 



no 



DO 70 

I = 1 ,NOP 


• 


m nrr-a— n r 

!22 a u - 


ID = 0 


112 

TTT 


13Y ♦ 30 








1 IM 
115 


DO 60 0 » 1 »K 

IFl IFIX(PPAKS(J»1 ) ) .NE* D 60 TO 60 



1 16 


CALL M0VAB5 10 » I 1 I) 



1 17 


CALL ANMODE 



118 


CALL DMPBUF 



1 19 


WRITE(6,55) alpha 1 oJ fPPARS 1 J t2 1 



120 

55 

FORPATl* » »6X , A6tF8»2 } 



121 


Y1 * 111 - 11. ' 



122 


CALL MOVEAC 10» »T1 ) 



. 123 


IF tlD .EQ. 0) GO TO 810 



128 


IDl » 0 



125 


102 c .0 



126 


call 0SHARC«J08.,Vl,ID»IDl«ID2|lDUP) 



1 27 


60 TO 811 



1 28 

810 

CALL DRAWA « 108 . t.T 1 ) 



129 

811 

CALL M0VAB5(0.I22) 



1 30 


CALL ANPODE 



1-3 1 


CALL DMPBUF 



132 


HRITE<6,55) ALPHA(O) iPPARS ( J,3) 



133 


Y1 = 122 ^ 13. 



1 38 


CALL MOVEaI lO.iTl ) 



135 


I F ( ID .EQ . 0 ) GO TO 8 13 



136 


10 1 =0 



137 


102 s 0 



138 


CALL 0SHARC{108.,YI,ID»ID1»ID2,IDUPJ 



139 


GO TO 815 



180 

813 

CALL DRAVIA ( 1 08, , Y 1 1 



191 

8 15 

10 * ID + 1 



142 


111 n 111 - 21 



183 


T22 a 122 ♦ 21 



188 

60 

CONTINUE 



1 85 


I I * I 1 - 13Y 




1M& 7D CONTINUE 

1*»7 X » (RUNSTP/ T TSPT ) + .OnOGtll 


1H8 I « ! 

C IF ( X .LT. I . ) I = 1 


150 

C 

IF < X ,GE, 60 . ) I ■ 3 

151 

C 

IFIX.6E.188U.JI«8 

152 


X o X V TSTEPA ( I ) 

153 


IF C5TARTX.LT-1 * 1 5 7 A R T X » T LO T B F ( I ) 

158 


X2 a ST ARTX / TSTEPX { I ) 

IBS 


X 1 » X2 ♦ X 

156 

157 


11 = TTSPt 
DO 75 J ■ 1 » 1 1 

T51 


XNON> = XI 

159 


X ! = X 1 + X 

1 *0 

75 

CONTINUE 

161 


CALL MOVABSIO, 100) 

1*2 


CALL ANMODE 

163 


CALL DMPBUF 



IFITTSPT ,6T. 8.) GO TO 80 

165 


irlTTSPT ,GT . .3. ) 60 TO 79 



WRITE(6»931 X2 ♦ 1 XNO { J )'i j» 1 *3 ) 

167 

93 

FORMAT!* •»15X»F6.1iUX,F6.I»H-X,F6.1,I0X,F6,1) 

TS1B 


GO TO 85 

169 

79 

1ARITE<6,77) X2i <XN01 J) »J»1 »8) 

170 

.77 

FORMAT!* *tl5X,F6.1,AX,F6,l,7X,F6.1,6X.F6,l,7X*F6,l) 

















86 FORWAT(» '.aOXt’tlML { • , A <4 » » ) • ) 

initialize dash INFnRRAT ion. 


DO 90 J » 1 >4 

I IDAS ( j» n = 0 


90 CONTINUE 

IFllEXECl »Eq. IPLL) RETURN 


VARIABLES* . • 
28 IL = 750 


IS2 a IS 
12 = IPLTPT 


00 90 I = 1 f NOP 
IL =» IL - I3T 


IDAS « -I 
IDUP * 1 


DO 30 n a I «K 

IF ( I F I X ( PPARS < II . 1) » *NE. J) GC TO 


lOAS = IDAS + 1 
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